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surface  testing.  The  baseline  smooth  wall  heat  transfer  data  show  that  the  qualifica¬ 
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Constant  wall  temperature  Stanton  numbers  determined  using  the  energy  balance 
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reported  up  to  x-Reynolds  numbers  of  3.5  million.  To  the  authors'  knowledge,  no 
smooth  wall  Stanton  number  data  have  been  reported  previously  in  the  3.5-10  million 
range  covered  by  the  THTTF  data.  Stanton  numbers  are  determined  in  the  THTTF  to 
uncertainties  of  2-5%,  depending  on  flow  conditions. 

The  initial  rough  wall  data  taken  using  the  THTTF  are  for  a  surface  roughened 
with  hemispherical  elements  spaced  four  base  diameters  apart.  The  Stanton  number 
and  skin  friction  coefficient  data  are  compared  with  predictions  from  the  discrete 
element  method  using  unmodified  the  roughness  models  developed  previous  to  this 
effort.  The  skin  friction  predictions  show  excellent  agreement  with  the  data.  This 
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is  encouraging  that  the  heat  transfer  predictions  are  as  good  as  they  are.  Further 
development  of  the  roughness  element  Nusselt  number  model  will  be  made  when  data 
are  taken  on  additional  rough  surfaces  in  the  THTTF. 
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SECTION  l 
INTRODUCTION 


Heat  transfer  and  skin  friction  can  be  significantly  larger 
for  a  turbulent  flow  over  a  rough  surface  compared  with  an  equiva¬ 
lent  turbulent  flow  over  a  smooth  surface.  Many  surfaces  of  engi¬ 
neering  interest  are  rough  in  the  aerodynamic  sense.  Missiles, 
re-entry  vehicles,  externally  mounted  stores  on  high  performance 
aircraft,  ship  hulls,  turbine  blades,  heat  exchangers,  piping 
networks  and  atmospheric  flows  are  examples  of  systems  in  which 
surface  roughness  can  play  an  important  role  in  heat  transfer  and 
skin  friction.  In  light  of  this  broad  applicability,  there  is 
significant  engineering  interest  in  the  development  of  accurate 
predictive  models  for  heat  transfer  and  fluid  mechanics  in  turbu¬ 
lent  flow  over  rough  surfaces. 

1 . 1  BACKGROUND 

Given  the  geometry  of  an  object  Immersed  In  a  flowfield,  a 
specification  of  the  freestream  flow  conditions,  and  a  geometrical 
description  of  the  roughness  of  the  system  surfaces,  an  analyst  or 
designer  would  like  at  least  to  be  able  to  predict  the  surface 
shear  distribution,  the  heat  transfer  distribution,  and  the  total 
drag.  In  the  past,  most  of  the  research  effort  was  to  develop 
computational  methods  for  various  geometries  with  smooth  surfaces, 
and  the  roughness  problem  has  received  relatively  little  attention. 
However,  many  systems  of  engineering  interest  have  surfaces  which 
are  aerodynamical ly  rough.  Therefore,  if  the  heat  transfer  and 
fluid  dynamics  are  to  be  predicted,  computational  procedures  to 
model  the  effects  of  rough  surfaces  must  be  developed  and  proven  by 
comparison  with  well -documented  data  sets. 

Schllchting  1936  [1]  experimentally  Investigated  the  fluid 
dynamics  of  this  type  of  problem.  He  related  his  skin  friction 
results  on  a  range  of  well-described  rough  surfaces  to  the  previous 
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results  obtained  oy  Nlkuradse  1933  [2]  for  sand-roughened  pipes 
through  definition  of  an  equivalent  sandgrain  roughness.  In  subse¬ 
quent  surface  roughness  effects  Investigations,  workers  used  these 
results  of  Schlichtlng  and  the  equivalent  sandgrain  roughness 
concept  to  analyze  their  experimental  data  and  to  develop  analyti¬ 
cal  models  for  use  in  predictive  methods. 

Over  the  past  decade  or  so,  a  predictive  approach  called  the 
discrete  element  method,  which  does  not  use  the  equivalent  sand¬ 
grain  roughness  concept,  has  been  used  with  varying  degrees  of 
rigor  by  several  groups  of  researchers  (Flnson  1975  [3],  Adams  and 
Hodge  1977  [4],  Finson  and  Wu  1979  [5],  Finson  and  Clark  1980  [6], 
Lin  and  Bywater  1980  [7],  Finson  1982  [8],  Coleman,  Hodge,  and 
Taylor  1983  [9],  Taylor,  Coleman  and  Hodge  1985  [10]).  Such 
approaches  rely  on  empirical  input  in  order  to  calibrate  the  rough¬ 
ness  models. 

Most  of  the  roughness-influenced  turbulence  data  taken  over 
the  years  has  been  on  ill-defined  rough  surfaces,  with  the  reported 
results  having  equivalent  sandgrain  roughness  values  implicitly 
included  at  some  stage  of  the  data  reduction.  Recently,  Coleman, 
Hodge,  and  Taylor  1983  [9]  found  that  Schllchting  had  made  errone¬ 
ous  assumptions  during  his  data  reduction  which  had  significant 
effects  on  the  data  which  he  reported.  The  reevaluation  of 
Schlichting's  data  (Coleman,  Hodge,  and  Taylor  1984  [11])  showed 
that  his  skin  friction  results  were  too  large  by  amounts  ranging  up 
to  73%  and  that  his  reported  equivalent  sandgrain  roughness  values 
were  high  by  amounts  ranging  from  26%  to  555%.  These  findings 
caused  some  consternation  since  practically  all  work  since  the 
1930's  on  surface  roughness  effects  relied  significantly  on  either 
the  skin  friction  or  equivalent  sandgrain  roughness  results  as 
originally  reported  by  Schllchting. 

A  series  of  data  sets  taken  at  Stanford  University  (Healzer 
1974  [12],  Plmenta  1975  [13],  Coleman  1976  [14],  Llgranl  1979  [15]) 
have  been  the  only  ones  reported  for  a  well-defined  rough  surface 
that  contain  heat  transfer  and  skin  friction  distributions  and 
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velocity,  temperature  and  Reynolds  stress  profiles.  However,  these 
data  sets  are  for  a  single  rough  surface  comprised  of  spheres  of  a 
single  size  packed  In  the  most  dense  array.  Holden  1983  [16]  re¬ 
ported  heat  transfer  and  skin  friction  distribution  measurements  on 
well-defined  surface  roughness  on  cones,  but  the  data  were  taken  at 
hypersonic  flow  conditions. 

Considering  the  limited  range  of  previously  reported  rough 
surface  data,  it  became  apparent  that  there  existed  a  critical  need 
for  accurate,  precise,  comprehensive  data  sets  on  both  the  heat 
transfer  and  the  fluid  dynamics  in  turbulent  flow  over  well-defined 
rough  surfaces.  It  was  concluded  that  if  a  reasonable  predictive 
capability  is  to  be  developed,  then  the  additional  experimental 
information  (particularly  for  heat  transfer)  must  be  obtained. 
Recognition  of  this  need  led  to  funding  under  the  DOD-Uni versity 
Research  Instrumentation  Program  of  a  water  tunnel  facility  and  the 
Turbulent  Heat  Transfer  Test  Facility  (THTTF)  in  the  Thermal  4 
Fluid  Dynamics  Laboratory  at  Mississippi  State  University.  A 
comprehensive  investigation  of  the  effects  of  surface  roughness 
element  size,  spacing  and  shape  on  skin  friction  in  fully  developed 
pipe  flow  over  a  wide  range  of  Reynolds  numbers  was  recently  com¬ 
pleted  using  the  water  tunnel  facility  in  this  laboratory  [17-21]. 

The  research  program  discussed  in  this  report  is  intended  to 
provide  accurate,  comprehensive  sets  of  measurements  of  Stanton 
number  distributions,  skin  friction  coefficient  distributions,  and 
profiles  of  velocity,  temperature  and  Reynolds  stresses  in  turbu¬ 
lent  boundary  layer  flows  over  surfaces  roughened  with  well-defined 
roughness  elements.  These  data  will  be  used  to  Improve  and  extend 
the  roughness  models  in  the  discrete  element  prediction  method, 
thus  expanding  our  capability  to  predict  the  effects  of  surface 
roughness  on  turbulent  flow  and  heat  transfer. 

The  current  phase  of  the  research  program  will  use  six  differ¬ 
ent  test  surfaces  in  the  THTTF— one  smooth  and  five  rough.  The 
smooth  surface  tests  serve  as  a  qualification  check  on  the  test 
facility  and  experimental  procedures.  Three  of  the  rough  test 
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surfaces  consist  of  smooth  plates  roughened  with  hemispherical 
roughness  elements  of  1.27  mm  diameter  spaced  2,  4,  and  10  base 
diameters  apart,  respectively.  The  other  two  rough  test  surfaces 
consist  of  smooth  plates  roughened  with  elements  that  are  truncated 
right  circular  cones  of  1.27  rrm  base  diameter  with  spacings  of  2 
and  4  base  diameters,  respectively. 

The  experimental  plan  is  to  test  each  set  of  surfaces  with 
basic  boundary  conditions  of  zero  pressure  gradient  and  constant 
wall  temperature  over  a  number  of  freestream  velocities  between 
about  6  and  67  m/s  such  that  the  total  set  of  data  will  thoroughly 
cover  behavior  in  the  aerodynamically  smooth,  transitionally  rough, 
and  fully  rough  regimes.  The  heat  transfer  data  are  the  most 
critical,  since  the  only  such  data  currently  available  on  a  well- 
described  rough  surface  are  those  from  the  previously  discussed 
Stanford  experiments,  which  used  a  single  rough  surface  composed  of 
1.27  mm  diameter  spheres  packed  in  the  most  dense  array. 

The  THTTF  was  designed  so  that  the  24  test  plates,  which  make 
up  its  test  surface,  can  be  replaced  with  a  new  set  without  com¬ 
pletely  tearing  down  the  test  section.  Some  re-instrumentation, 
calibration  and  qualification  is  necessary  for  each  new  set  of  test 
plates,  however,  to  maintain  the  high  accuracy  and  precision  which 
are  an  inherent  part  of  the  overall  objective  of  this  test  program. 

1.2  OVERVIEW  OF  REPORT 

This  report  presents  results  for  the  smooth  test  surface  and  a 
rough  surface  with  hemispherical  elements  at  a  spacing  of  four  base 
diameters.  The  first  tests  were  the  calibration,  qualification  and 
general  " shake-down "  of  the  facility  using  the  smooth  test  surface. 
The  ability  to  reproduce  accepted  smooth  wall  results  for  non- 
isothermal  turbulent  boundary  layer  flow  in  the  THTTF  was  deemed 
necessary  before  proceeding  with  rough  wall  investigations  using 
the  facility.  In  addition,  the  smooth  wall  data  provide  an  appro¬ 
priate  baseline  with  which  to  compare  the  data  obtained  for  rough 
walls  using  the  same  test  apparatus. 
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This  report  Is  organized  In  the  following  manner.  In  Section 
2,  the  THTTF  and  each  of  Its  four  sub-systems  are  described.  In 
Section  3,  smooth  wall  fluid  dynamics  data  Including  skin  friction 
coefficient  distributions  and  boundary  layer  profiles  of  mean 
velocity  and  Reynolds  stresses  are  presented  and  compared  to  ac¬ 
cepted,  previously  reported  smooth  wall  data.  The  smooth  wall 
qualification  heat  transfer  data,  Including  Stanton  number  distri¬ 
butions  and  boundary  layer  temperature  profiles,  are  presented  and 
compared  to  accepted,  previously  reported  smooth  wall  data  In 
Section  4.  Also  presented  are  Stanton  number  data  which  extend  the 
previously  published  smooth  wall  data  from  x-Reynolds  numbers  of 
3.5  million  to  10  million.  In  Section  5,  additional  smooth  wall 
heat  transfer  results  are  presented  for  a  variety  of  non-constant 
wall  temperature  boundary  conditions.  In  Section  6,  the  fluid 
dynamics  and  heat  transfer  data  from  the  rough  wall  are  presented 
and  compared  with  calculations  from  the  discrete  element  prediction 
method.  The  results  from  the  research  program  are  summarized  In 
Section  7. 

Detailed  discussions  of  the  experimental  determination  of 
Stanton  number,  Its  uncertainties,  and  estimates  of  the  uncertain¬ 
ties  in  the  various  experimental  measurements  are  presented  in 
Appendix  I.  Probe  measurements  and  calibrations  are  described  In 
Appendix  II  for  the  hot-wire  and  thermocouple  probes,  while  the 
determination  of  skin  friction  coefficients  Is  covered  in  Appendix 
III. 
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SECTION  2 

DESCRIPTION  OF  EXPERIMENTAL  FACILITY 


The  Turbulent  Heat  Transfer  Test  Facility  (THTTF)  was  con¬ 
structed  for  the  experimental  Investigation  of  heat  transfer  and 
fluid  dynamics  behavior  in  turbulent  boundary  layers  over  rough 
surfaces.  The  design  of  the  THTTF  was  based  on  the  preliminary 
analysis  and  design  work  of  Norton  1983  [22].  The  THTTF  is  geometri¬ 
cally  similar  to  the  test  apparatus  used  in  the  Stanford  University 
program  that  investigated  turbulent  boundary  layer  flow  and  heat 
transfer  in  flow  over  a  single,  porous  rough  surface  with  transpi¬ 
ration  (Healzer  1974  [12],  Pimenta  1975  [13],  Coleman  1976  [14], 
Ligrani  1979  [15]). 

2.1  GENERAL  DESCRIPTION 

The  THTTF,  shown  schematically  in  Figure  1,  is  a  closed  loop 
subsonic  wind  tunnel  designed  to  deliver  a  uniform  air  flow  over  a 
set  of  24  individually  heated  flat  test  plates  which  are  abutted 
together  to  form  a  continuous  flat  surface.  Each  plate  can  be 
maintained  at  a  constant  uniform  temperature,  and  each  set  of  rough 
plates  has  a  well  defined  surface  roughness.  The  Stanton  number 
distribution  along  the  test  surface  is  obtained  by  applying  an 
energy  balance  to  each  test  plate.  Distribution  of  the  local  skin 
friction  coefficient  along  the  test  surface  and  other  boundary 
layer  parameters  are  determined  with  hot-wire  anemometry.  Thermal 
boundary  layer  parameters  are  determined  from  temperature  profiles 
measured  with  a  thermocouple  probe  and  the  appropriate  hot-wire 
data. 

Measurement  of  individual  variables  such  as  air  velocity  and 
plate  temperature,  the  calibration  of  the  Instruments  used  to 
measure  these  variables,  the  determination  of  experimental  Stanton 
numbers  from  these  measured  variables,  and  an  analysis  of  the 
uncertainty  associated  with  these  Stanton  numbers  are  presented  in 
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Appendix  I.  The  boundary  layer  probe  measurement  techniques  are 
giver  In  Appendix  II,  and  the  methods  of  determining  local  skin 
friction  coefficients  are  discussed  in  Appendix  III. 

Four  primary  systems  are  required  for  the  proper  operation  of 
the  THTTF :  (1)  the  air  flow  system,  (2)  the  test  plate  system,  (3) 
the  cooling  water  system,  and  (4)  the  Automated  Data  Acquisition 
and  Control  System  (ADACS).  These  are  discussed  in  detail  below. 

2.2  AIR  FLOW  SYSTEM 

The  air  flow  system  is  a  closed  loop  system  designed  to  de¬ 
liver  a  uniform  velocity  (6  to  67  m/s),  low  turbulence  intensity, 
controlled  temperature  air  flow  at  the  10.2  by  50.8  cm  (4.0  by  20.0 
in)  inlet  of  the  2.4  m  (8.0  ft)  long  test  section  which  contains 
the  test  surface.  These  air  velocities  correspond  to  an  x-Reynolds 
number  range  of  about  one  million  to  ten  million  at  the  downstream 
end  of  the  test  section.  A  Buffalo  Forge  size  45AW  industrial 
blower  is  the  prime  mover  for  the  air  flow  system.  The  blower  has 
a  rating  of  198  cubic  meters  of  air  per  minute  at  38  cm  of  water 
static  pressure.  An  18.6  kilowatt  (25-hp)  Dynamatic  Model  ACM-280 
electric  motor  with  an  Eaton  variable  speed  eddy  current  clutch 
drive  system  drives  the  blower  with  a  belt  and  pulley  system.  Air 
velocity  in  the  test  section  is  set  using  the  Eaton  eddy  clutch 
controller  to  control  the  rotation  speed  of  the  blower.  The  con¬ 
troller  can  be  adjusted  and  set  using  a  manually  adjusted  potenti¬ 
ometer  or  a  dc-voltage  signal  from  the  ADACS. 

Air  exiting  the  blower  enters  a  1.2  m  wide  by  0.6  m  tall 
wooden  overhead  duct  which  is  connected  to  the  blower  and  header  by 
flexible  couplings.  The  air  then  turns  through  the  header  and 
passes  through  a  linen  cloth  filter  in  the  filter  box.  Next  the 
air  passes  through  an  air/water  heat  exchanger  with  a  4  row  cooling 
coll.  Upon  leaving  the  heat  exchanger,  the  air  passes  through  a  3.8 
cm  thick  aluminum  honeycomb  with  a  cell  length-to-diameter  ratio  of 
6  and  then  through  a  series  of  4  woven  stainless  steel  screens  with 
an  open  area-ratio  of  0.598  and  a  wire  diameter  of  0.136  mm. 
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Following  the  screens,  the  air  enters  a  three-dimensional, 

19.8  to  1  contraction  ratio  fiberglass  nozzle  with  a  84  cm  by  122 
cm  (33  In  by  48  In)  Inlet  and  a  10  cm  by  51  cm  (4  In  by  20  In) 
outlet.  The  nozzle  was  designed  (Healzer  1974  [12])  to  smoothly 
accelerate  the  flow  without  separation  at  the  nozzle  Inlet  or 
outlet.  Uniform  velocity  air  is  delivered  from  the  nozzle  to  the 
test  section  inlet.  Measurements  at  freestream  air  velocities  of 
12  and  28  m/s  indicated  the  axial  velocity  at  the  nozzle  exit  is 
uniform  within  about  0.5%.  Freestream  turbulence  intensities 
measured  4  cm  downstream  of  the  nozzle  exit  were  less  than  0.3%  for 
freestream  velocities  from  61  m/s  down  to  6  m/s  and  less  than  0.4% 
at  3  m/s. 

The  test  section  has  clear  cast  acrylic  (plexiglass)  sidewalls 
of  1.3  cm  thickness  and  an  adjustable,  flexible,  clear  cast  acrylic 
top  of  1.3  cm  thickness.  The  bottom  wall  is  made  up  of  the  test 
plates.  Static  pressure  taps,  spaced  every  10.2  cm  even  with  the 
center  of  each  plate  in  the  test  surface,  are  located  in  one  of  the 
sidewalls  approximately  2.5  cm  above  the  plates.  Access  holes  for 
test  probes  are  located  along  the  center  of  the  top  wall,  centered 
over  each  plate,  and  transversely  at  strategic  locations.  The 
holes  are  plugged  with  precision  machined  acrylic  stoppers  when 
probe  access  is  not  required. 

The  flexible  upper  surface  of  the  test  section  can  be  adjusted 
to  maintain  the  prescribed  zero  pressure  gradient  along  the  flow 
direction.  A  Dwyer  inclined  water  manometer  with  a  resolution  of 
0.06  mm  of  water  was  connected  to  the  static  pressure  taps  along 
the  sidewall  of  the  test  section  to  measure  the  pressure  gradient 
In  the  flow  direction.  The  difference  in  static  pressure  In  the 
test  section  between  all  pressure  tap  locations  and  that  tap  lo¬ 
cated  at  the  second  plate  was  maintained  at  less  than  0.19  mm  of 
water  for  a  velocity  of  12  m/s  and  to  less  than  0.32  mm  of  water 
for  a  velocity  of  43  m/s. 
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As  it  exits  the  test  section,  the  air  passes  through  an  ad¬ 
justable  plexiglass  diffuser  which  links  the  test  section  to  a 
vaned  wooden  diffuser.  Easily  removed  screen  Inserts,  which  are 
used  to  produce  a  pressure  drop  and  thus  used  as  a  coarse  adjust¬ 
ment  of  the  absolute  static  pressure  level  In  the  test  section,  are 
located  at  the  entrance  of  the  wooden  diffuser.  Following  the 
vaned  diffuser,  the  air  enters  the  blower  plenum  from  which  it 
passes  through  a  flexible  coupling  Into  the  blower  intake. 

Suction  and  ejection  of  air  to  and  from  the  test  section 
through  any  small  air  gaps  In  the  test  section  were  minimized  by 
equalization  of  the  test  section  static  pressure  with  the  ambient 
pressure.  Filtered  make-up  air  to  replace  the  inevitable  air 
leakage  from  the  overhead  ducting,  filter  boxes,  and  heat  exchanger 
was  ducted  through  a  box  of  adjustable  orifices  to  the  blower 
plenum.  Very  fine  balancing  of  the  test  section  air  static  pres¬ 
sure  and  the  ambient  pressure,  typically  to  within  ±0.13  mm  of 
water,  was  accomplished  by  adjusting  the  orifices  in  the  make-up 
air  box. 

In  order  to  minimize  vibrations,  the  blower  and  blower 
motor,  which  are  the  primary  sources  of  mechanical  vibrations, 
are  mounted  on  a  massive  concrete  pad  with  vibration  dampening 
feet.  Transfer  of  mechanical  vibrations  throughout  the  THTTF  has 
been  minimized  by  the  use  of  non-rigid  joints  at  key  locations  in 
the  flow  path  to  effectively  Isolate  the  test  section  from  vibra¬ 
tion  sources.  Noise  from  air  flowing  in  the  overhead  duct,  plenum, 
and  header  has  been  reduced  by  lining  these  air  passages  with  batt 
Insulation  covered  by  rigid  fiberglass  insulation  board.  Blower 
and  blower  motor  noise  have  been  reduced  by  housing  the  blower  and 
motor  in  vented,  insulation-lined  boxes. 

2.3  PLATE  SYSTEM 

A  cross  sectional  view  of  the  test  section  Is  shown  In  Figure 
2.  The  test  plates  are  supported  on  precision  straight  edges  which 
are  thermally  Isolated  from  the  steel  side  rails  which  provide  the 
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primary  structural  support.  These  side  rails  are  heated  and  act  as 
guard  heaters  that  help  to  minimize  the  conduction  heat  losses  from 
the  plates.  A  removable.  Insulated,  wooden  enclosure  around  the 
base  of  the  test  section  reduces  heat  loss  from  the  metal  support 
rails  to  the  laboratory  environment. 

The  24  plates  which  comprise  the  test  surface  are  each  10.2  cm 
(4.0  In)  in  the  flow  direction  by  45.7  cm  (18.0  in)  in  the  trans¬ 
verse  direction  by  0.95  cm  (0.375  in)  thick.  The  precision  ma¬ 
chined  test  plates  are  made  of  electroless  nickel  plated  aluminum. 
The  smooth  surface  plates  used  in  the  baseline  tests  have  a  surface 
finish  with  centerline  average  roughness,  Ra,  measured  as  less  than 
0.5  micrometers.  The  rough  wall  plates  discussed  In  this  report 
have  a  centerline  average  roughness  measured  as  less  than  1.6 
micrometers  on  the  "smooth"  wall  portion  of  the  plates.  The  plates 
are  assembled  using  dowels  to  form  a  continuous  and  smooth  test 
surface  as  shown  in  Figure  3.  The  allowable  step  (or  mismatch)  at 
the  joint  between  two  plates  Is  0.013  cm. 

Each  plate  is  Instrumented  with  two  thermistors  for  tempera¬ 
ture  measurement,  and  each  has  its  own  motor-driven  variable  volt¬ 
age  transformer/plate  heater  circuit  which  is  controlled  by  the 
AOACS.  A  flexible  resistance  heating  pad  (plate  heater)  affixed  to 
the  bottom  of  each  plate  provides  a  uniform  heat  flux  to  the  lower 
plate  surface.  The  plate  heaters,  which  were  custom  manufactured  by 
Watlow  Electric  Manufacturing  Co.,  are  about  one  mm  thick  and  are 
made  of  resistance  wire  spiraled  around  a  glass  cord  sandwiched 
between  two  pieces  of  glass  fabric  coated  with  silicone  rubber. 

All  heater  circuits  are  supplied  electrical  power  through  the 
same  Powermark-75110  A.C.  voltage  regulator  which  is  connected  to 
the  building  service  (110  volt)  as  shown  In  Figure  4.  Fine 
adjustment  of  the  plate  heater  power  Is  accomplished  with  a 
Powerstat-15M21  motor  driven  variable  transformer  In  each  heater 
circuit.  These  motor  driven  transformers  are  grouped  In  banks  of 
8,  which  are  supplied  power  through  1  of  3  manually  set  Varlac-WlO 
variable  transformers  used  for  gross  step-down  of  the  regulated 
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power.  Experience  In  acquiring  the  heat  transfer  data  has  shown 
that  plate  temperatures  can  be  held  within  ±0.1  C  of  a  prescribed 
constant  temperature  boundary  condition  along  the  entire  test 
section. 

2.4  COOLING  WATER  SYSTEM 

The  cooling  water  system  maintains  the  test  air  at  a  constant 
temperature  by  extracting  energy  added  by  the  heated  plates  and 
blower.  Cool  water  from  a  568  liter  water  storage  tank  Is  moved  by 
a  186  watt  (1/4  hp)  Bell  &  Gossett  Model  1522  pump  through  5  cm  pvc 
piping  to  a  Trane  air/water  heat  exchanger.  The  heat  exchanger  has 
4  rows  (in  the  flow  direction)  of  finned  cooling  coils  with  a  84  cm 
by  112  cm  (33  In  by  48  In)  area  normal  to  the  flow.  Heat  from  the 
test  air  is  transferred  to  the  cooling  water  via  the  cooling  coils 
before  the  air  enters  the  nozzle  and  then  the  test  section.  Pro¬ 
portional  amounts  of  warmer  water  returning  from  the  cooling  coils 
to  the  storage  tank  can  be  dumped  into  a  floor  drain  through  a  1.9 
cm  motorized  ball  valve,  GF-Type  105.  This  ball  valve  may  be 
adjusted  and  controlled  manually  or  by  the  ADACS.  Water  level  in 
the  storage  tank  Is  maintained  at  a  desired  level  with  make-up 
cooling  water  dispensed  through  adjustable  depth  sensing  valves  fed 
by  the  building  supply. 

2.5  DATA  ACQUISITION  SYSTEM 

A  Hewlett-Packard  Series  9000  Model -220  microcomputer  and  a 
Hewlett-Packard  3054A  Automatic  Data  Acquisition  and  Control  System 
are  used  to  monitor  and  control  the  THTTF.  The  ADACS  Includes  an 
HP-3437A  high  speed  system  voltmeter,  an  HP-3456A  high  resolution 
digital  voltmeter,  an  HP-3497A  data  acquisition/control  unit  and  a 
number  of  special  function  plug-in  assemblies.  A  detailed  discus¬ 
sion  of  the  ADACS  and  Its  use  in  the  THTTF  Is  given  by 
Suryanarayana  1986  [23]. 
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Transducers  that  monitor  the  THTTF  are  wired  Into  the  ADACS, 
which  relays  the  Information  to  the  microcomputer.  The  microcompu¬ 
ter  digests  the  operating  condition  Information,  decides  on  the 
proper  response  based  on  programming,  and  sends  commands  for  the 
proper  controller  response  to  the  ADACS.  The  ADACS  can  control  the 
rail  heaters,  plate  heaters,  cooling  system  dump  valve,  and  blower 
motor  in  response  to  instructions  from  the  microcomputer.  When  the 
THTTF  has  been  brought  to  the  desired  equilibrium  conditions,  the 
computer  directs  the  ADACS  to  perform  the  necessary  data  collec¬ 
tion. 


12 


Figure  1.  Schematic  of  the  Turbulent  Heat  Transfer  Test  Facility  (THTTF) . 
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Approximately  to  scale. 


Figure  4.  Schematic  of  the  Plate  Heater  Power  Circuit 
and  Control  Loop. 
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SECTION  3 

SMOOTH  WALL  QUALIFICATION:  FLUID  DYNAMICS 


The  Turbulent  Heat  Transfer  Test  Facility  (THTTF)  was  designed 
to  provide  quality  heat  transfer  data  primarily  for  zero  pressure 
gradient,  constant  wall  temperature,  incompressible  flow  over  flat 
plates  with  various  surface  roughnesses.  To  produce  quality  con¬ 
vective  heat  transfer  data,  the  THTTF  must  produce  flows  with  prop¬ 
erties  within  acceptable  fluid  dynamics  bounds,  and  the  techniques 
used  in  collection  of  the  fluid  dynamics  data  must  be  proven.  The 
discussion  in  this  section  is  intended  to  show  that  the  THTTF  is 
operating  within  acceptable  fluid  dynamics  bounds  and  to  document 
the  validity  of  the  Instrumentation,  data  collection  and  data 
reduction  procedures.  This  objective  is  met  by  comparing  the  smooth 
wall  data  obtained  in  the  THTTF  with  previously  published,  well- 
accepted  smooth  wall  results  of  others.  All  comparisons  in  this 
section  are  for  zero  pressure  gradient,  incompressible,  isothermal 
turbulent  boundary  layer  flow  over  a  smooth  surface. 

Fluid  dynamics  variables  Including  profiles  of  mean  velocity 
and  Reynolds  stresses  and  distributions  of  local  skin  friction 
coefficients  were  determined.  A  horizontal  hot-wire  probe  and  a 
slanted  hot-wire  probe  were  used  to  measure  mean  velocity  and 
turbulence  quantities  as  discussed  in  Appendix  II.  As  discussed  in 
Appendix  III,  skin  friction  coefficients  were  determined  from 
Reynolds  shear  stress  measurements  made  using  hot-wire  anemometry 
and  also  from  Preston  tube  measurements.  In  the  following,  data 
taken  using  the  THTTF  are  compared  with  data  obtained  on  other  test 
facilities  which  have  generated  definitive,  accepted  data. 

3.1  FREESTREAM  FLOW  QUALITY 

Pitot  tube  surveys  of  the  air  flow  entering  the  test  section 
were  made  at  the  nozzle  exit  for  freestream  air  velocities  of  12 
and  28  m/s.  The  results  Indicate  that  the  flow  Is  uniform  to  within 
about  0.5%  of  the  mean  velocity. 
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As  described  previously,  screens  and  honeycomb  are  used  Imme¬ 
diately  downstream  of  the  heat  exchanger  to  suppress  turbulent 
fluctuations  In  the  air  flow  entering  the  nozzle.  Hot-wire 
anemometrv  was  used  to  measure  the  freestream  turbulence  intensity 
with  the  horizontal  wire  probe  positioned  at  the  test  section 
centerline  4  cm  downstream  of  the  nozzle  exit.  These  measurements 
showed  the  freestream  turbulence  intensity  to  be  less  than  0.3%  for 
freestream  velocities  from  61  m/s  down  to  6  m/s  and  less  than  0.4% 
at  3  m/s.  This  compares  favorably  with  freestream  turbulence 
intensities  between  2  and  5  percent  reported  by  Reynolds,  Kays  and 
Kline  1958  [24],  whose  smooth  wall  Stanton  number  data  are  consid¬ 
ered  definitive,  for  freestream  velocities  ranging  from  14  to  39 
m/s.  Freestream  turbulence  Intensity  of  about  0.7%  was  reported  by 
Moffat  1967  [25],  Kearney  1970  [26],  and  Orlando  1974  [27]  for  a 
wind  tunnel  somewhat  similar  In  design  to  the  THTTF,  and  Pimenta 
1975  [13]  reported  a  value  of  about  0.4%  for  the  Stanford  Roughness 
Rig  which  is  geometrically  similar  to  the  THTTF. 

3.2  BOUNDARY  LAYER  TWO  DIMENSIONALITY 

The  two  dimensionality  of  the  flow  was  investigated  by  taking 
velocity  profiles  with  the  horizontal  hot-wire  at  transverse  sta¬ 
tions  spaced  5  cm  apart  over  plate  11  (midway  down  the  length  of 
the  test  section  at  x  »  1.07  m).  Boundary  layer  velocity  profile 
sets  were  taken  at  freestream  velocities  of  12  and  27  m/s.  Momen¬ 
tum  thicknesses  determined  from  these  transversely  positioned 
velocity  profiles  are  shown  In  Figure  5.  The  variation  in  the 
momentum  thicknesses  across  plate  11  in  the  transverse  direction 
was  less  than  ±5%  about  the  mean  for  both  velocities.  Healzer  1974 
[12]  reported  a  variation  of  the  momentum  thickness  in  the  trans¬ 
verse  direction  of  less  than  3%  at  a  velocity  of  27  m/s  for  the 
Stanford  Roughness  Rig. 

An  additional  check  on  the  two  dimensionality  of  the  boundary 
layer  was  made  by  measuring  the  spanwlse  variation  In  the  skin 
friction  coefficient  at  plate  11  for  a  flow  velocity  of  12  m/s. 

Skin  friction  coefficients  determined  from  a  Preston  tube  measuring 
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at  positions  spaced  0.6  cm  (1/4  In)  transversely  across  plate  11 
varied  by  less  than  9%  peak  to  peak  and  less  than  5%  from  the  mean. 
Mehta  and  Hoffman  1987  [28]  reported  spanwise  variations  in  local 
friction  coefficients  for  their  tests  in  two  wind  tunnels  at  Impe¬ 
rial  College  between  8%  and  20%  peak  to  peak. 

3.3  SKIN  FRICTION  COEFFICIENT  RESULTS 

The  definitive  skin  friction  coefficient  data  for  zero  pres¬ 
sure  gradient  turbulent  boundary  layer  flows  over  smooth  surfaces 
are  those  of  Schultz-Grunow  1941  [29].  Schultz-Grunow  determined 
the  drag  across  flat  plates  in  turbulent  flows  of  near  zero  pres¬ 
sure  gradient  with  a  "direct"  measurement  method,  which  used  a 
sensitive  lever-arm  and  wire  torsion  system  to  measure  the  drag  on 
a  movable  plate  section.  He  curvefit  his  data  to  get  the  explicit 
expression 


Cf/Z  -  0.185(log10Rex)“2*584  (1) 

Schultz-Grunow' s  expression  has  been  shown  to  match  the  data  of 
numerous  experimenters  to  within  about  ±10%  and  is  referenced  as 
definitive  by  Schlichting  1979  [30]  and  Kays  and  Crawford  1980 
[31]. 

Equation  (1)  is  written  in  terms  of  the  distance  from  the 
plate's  leading  edge  (x)  for  an  Ideal  turbulent  boundary  layer 
starting  exactly  at  the  leading  edge.  A  flat  plate  boundary  layer 
that  is  tripped  turbulent  acts  as  If  the  turbulent  boundary  layer 
starts  at  a  fictitious  point  (virtual  origin)  upstream  of  the  trip. 
This  virtual  origin  corresponds  to  the  point  where  the  turbulent 
boundary  layer  would  begin  if  it  was  not  preceded  by  a  laminar 
boundary  layer.  Therefore,  experimental  skin  friction  data  from 
the  THTTF,  which  has  a  tripped  boundary  layer,  is  presented  on  the 
basis  of  a  virtual  origin  when  comparing  It  to  the  skin  friction 
correlation  in  terms  of  x. 
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The  expression  used  to  determine  the  virtual  origin  from 
measured  velocity  profile  data  at  each  profile  station  Is  that 
given  by  Kays  and  Crawford  as 


x|  -  0-036  (Rexv)“0-2 


(2) 


where  xv  is  the  streamwise  coordinate  measured  from  the  virtual 
origin.  For  a  given  freestream  velocity,  the  "best"  virtual  origin 
was  obtained  by  application  of  a  least  squares  approach  to  the 
differences  between  the  momentum  thicknesses  determined  from  the 
measured  sets  of  velocity  profiles  and  the  corresponding  momentum 
thicknesses  predicted  by  equation  (2)  for  assumed  virtual  origins. 
The  virtual  origins  determined  in  this  manner  were  5,  14,  and  27  cm 
upstream  of  the  nozzle  exit  for  freestream  velocities  of  12,  27, 
and  43  m/s,  respectively.  Values  of  35  and  43  cm  were  determined 
by  extrapolation  for  freestream  velocities  of  58  and  67  m/s. 

Measurements  of  skin  friction  coefficient  distributions  along 
the  test  surface  were  made  for  U„  »  12  and  43  m/s  using  the 
hot-wire  method  and  U_  *  12,  27,  43,  58  and  67  m/s  using  the 
Preston  tube  method.  Both  these  techniques  are  discussed  in  Appen¬ 
dix  III. 

Figure  6  shows  the  local  skin  friction  coefficient  distribu¬ 
tions  in  the  THTTF  determined  by  hot-wire  anemometry  compared  with 
the  Schultz-Grunow  expression  with  a  ±10%  range  Indicated.  The 
scatter  in  the  skin  friction  coefficient  data  determined  with 
hot-wire  anemometry  fall  essentially  within  the  ±10%  bands  about 
the  accepted  correlation.  This  is  especially  encouraging,  since 
this  is  the  sole  technique  used  In  this  research  program  to  deter¬ 
mine  skin  friction  coefficients  In  rough  surface  flows.  It  also 
serves  as  a  verification  check  on  the  estimated  10  to  12%  uncer¬ 
tainty  in  Cf  from  the  hot-wire  method. 

Skin  friction  coefficients  determined  by  the  Preston  tube 
method  are  shown  In  Figure  7  compared  with  equation  (1)  and  the 
±10%  bands.  A  composite  plot  of  the  skin  friction  data  determined 
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by  both  techniques  Is  given  In  Figure  8.  This  figure  shows  that 
essentially  all  of  the  present  data  agrees  with  the  well  accepted 
Schultz-Grunow  expression  within  the  ±10%  band.  Therefore,  It  can 
be  concluded  that  the  local  skin  friction  coefficient  data  from  the 
THTTF  is  in  substantial  agreement  with  the  definitive  skin  friction 
correlation  and  that  the  facility  and  measurement  techniques  are 
performing  correctly  for  skin  friction  measurements. 

The  momentum  thickness  (62)  of  a  boundary  layer  accounts  for 
the  upstream  history  of  the  flow  to  a  considerable  degree.  Thus, 
skin  friction  data  from  the  THTTF  are  also  compared  with  an  ac¬ 
cepted  skin  friction  correlation  based  on  the  S2-Reynolds  number. 
For  zero  pressure  gradient,  incompressible  flow  over  smooth  flat 
plates,  Kays  and  Crawford  recommend 


Cf/2  -  0.0125{Refi2)"0*25  (3) 

Figure  9  shows  the  local  skin  friction  coefficient  distributions 
determined  by  hot-wire  anemometry  versus  the  52-Reynolds  number 
computed  from  the  measured  momentum  thickness  at  each  profile 
station  and  also  presents  equation  (3)  with  a  ±10%  range  indicated. 
Skin  friction  coefficients  determined  by  the  Preston  tube  method 
are  shown  versus  the  52-Reynolds  number  in  Figure  10  along  with 
equation  (3)  and  the  ±10%  bands.  Fewer  Preston  tube  determined 
skin  friction  coefficients  are  shown  in  Figure  10  than  in  Figure  7 
because  Preston  tube  measurements  were  made  at  more  locations  than 
were  velocity  profiles  from  which  momentum  thicknesses  were  deter¬ 
mined. 

A  composite  plot  of  the  skin  friction  data  as  determined  by 
both  methods  is  given  in  Figure  11.  Inspection  of  this  figure 
shows  that  essentially  all  of  the  data  agrees  with  equation  (3) 
within  a  ±10%  band.  This  comparison  provides  additional  confidence 
in  the  skin  friction  measurements. 


21 


3.4  BOUNDARY  LAYER  PROFILE  DATA 

Boundary  layer  profiles  of  mean  velocity  and  turbulence  quan¬ 
tities  were  also  measured  for  comparison  with  "standard"  smooth 
wall  results.  Boundary  layer  mean  velocity  profiles  taken  with  the 
horizontal  hot-wire  over  test  plate  19  (x  -  1.88  m)  for  freestream 
velocities  of  12  and  43  m/s  are  shown  in  Figure  12  plotted  in  u/U. 
vs  y/S  coordinates.  These  profiles  are  plotted  In  Figure  13  In 
inner  region  coordinates  (u+  vs  y+)  and  are  compared  with  the  usual 
smooth  surface  law  of  the  wall  expression 

u+  -  2.44  in  y+  +  5.0  (4) 

The  friction  velocity  used  In  u+  and  y+  was  that  determined  by  the 
hot-wire  method. 

Profiles  of  the  axial  direction  turbulence  fluctuations  were 
measured  with  a  horizontal  hot-wire  at  6  stations  spaced  40.6  cm 
apart  along  the  length  of  the  test  section  as  discussed  In  Appendix 
II.  These  profiles  at  plates  15,  19,  and  23  for  U„  -  12  m/s  are 
shown  in  Figure  14  In  Inner  region  coordinates  and  are  compared 
with  curves  representing  the  axial  turbulence  data  of  Klebanoff 
1955  [32]  and  Laufer  1954  [33].  Figure  15  shows  the  corresponding 
profiles  for  U.  -  43  m/s.  The  data  shown  In  these  two  figures 
collapse  fairly  well  In  the  Inner  region  and  follow  the  trend  of 
the  data  of  klebanoff  and  Laufer. 

These  same  profiles  for  U„  -  12  m/s  are  shown  versus  y/S  In 
Figure  16.  The  data  follow  the  trend  of  Klebanoff' s  data,  repre¬ 
sented  by  the  solid  curve.  In  the  outer  portion  of  the  boundary 
layer,  the  axial  turbulence  fluctuations  In  the  THTTF  are  greater 
than  those  of  Klebanoff  since  his  tunnel  had  a  lower  freestream 
turbulence  intensity  (0.02%  for  U.  -  9  m/s  and  0.04%  for  U.  -  31 
m/s).  Corresponding  comparisons  are  shown  In  Figure  17  for  U„  -  43 
m/s.  The  nondlmenslonal  axial  turbulence  distributions  In  the 
THTTF  at  12  and  43  m/s  are  almost  Identical  except  In  the  lower  5% 
of  the  boundary  layer.  The  horizontal  hot-wire  length  was 
greater  than  100  wall  units  (v/u*)  at  the  higher  freestream 
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velocities.  Thus,  as  discussed  by  Llgranl  and  Bradshaw  1987  [34], 
the  hot-wire  used  In  this  study  was  unable  to  accurately  resolve 
the  axial  turbulence  very  near  the  surface  at  higher  freestream 
velocities. 

__The  profiles  of  all  three  turbulence  fluctuating  quantities 
((ui2)l/2/u*,  (v ' 1/2/u* ,  (w • 2) 1/2/u*)  were  measured  at  plate  11 
for  U„  -  12  m/s  by  traversing  the  boundary  layer  with  a  horizontal 
hot-wire  probe  and  a  45°  slant  hot-wire  probe  as  described  in 
Appendix  II.  These  data  are  shown  in  Figure  18  and  compared  with 
solid  curves  representing  the  data  of  Klebanoff.  The  profile  of 
the  Reynolds  shear  stress  at  the  same  position  and  conditions  Is 
shown  in  Figure  19  and  compared  with  a  solid  curve  which  represents 
the  Reynolds  shear  stress  data  of  Klebanoff.  The  Reynolds  shear 
stress  data  taken  on  the  THTTF  are  in  general  agreement  with  the 
curve  representing  Klebanoff1 s  data. 
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Figure  5.  Momentum  Thickness  Distributions  Measured  Transve 
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Grunow  (1941). 
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Figure  12.  Comparison  of  Mean  Velocity  Profiles  at  Plate  19  for  Nominal  Freestream 
Velocities  of  12  and  43  m/s. 
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Plate  15 
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Figure  16.  Axial  Turbulence  Intensity  Profiles  vs.  y/6  for  a  Nominal  Freestream  Velocity 
of  12  m/s  Compared  with  the  Data  of  Klebanoff  (1955). 


Plate  15,  x  =  1.47 


12  m/s 


Freestream  Velocity  of  12  m/s  Compared  with  the  Data  of  Klebanoff 


12  m/s 


Reynolds  Shear  Stress  at  Plate  19  for  a  Nominal  Freestream 
Compared  with  the  Data  of  Klebanoff  (1955). 


SECTION  4 

SMOOTH  WALL  QUALIFICATION:  HEAT  TRANSFER 


The  primary  objective  of  this  section  Is  to  present  the  smooth 
wall  Stanton  number  measurements  from  the  THTTF,  compare  them  with 
previously  published  data  from  accepted  sources,  and  verify  the 
fitness  of  the  test  rig  and  the  correctness  of  the  instrumentation, 
data  acquisition  system  and  data  reduction  procedures  used  In 
obtaining  the  Stanton  number  data  and  the  estimates  of  its  uncer¬ 
tainty.  All  comparisons  and  data  presented  in  this  section  are  for 
zero  pressure  gradient,  constant  wall  temperature,  incompressible 
boundary  layer  flow  over  smooth  flat  plates. 

4.1  STANTON  NUMBER  RESULTS 

Stanton  numbers  were  determined  from  energy  balances  on  each 
test  plate,  as  described  In  detail  in  Appendix  I.  The  THTTF  data 
are  presented  using  a  definition  of  the  Stanton  number  which  is 
based  on  the  difference  between  the  wall  temperature  and  the  free- 
stream  total  temperature.  Other  data  with  which  comparisons  are 
made  are  based  on  a  definition  of  the  Stanton  number  which  uses  the 
difference  between  the  wall  temperature  and  the  freestream  recovery 
temperature.  This  difference  In  Stanton  number  definitions  is  neg¬ 
ligible  in  the  comparisons  made  because  the  differences  in  the 
total  and  recovery  temperatures  are  numerically  insignificant  for 
the  range  of  air  velocities  considered  by  the  previous  experiment¬ 
ers. 

The  definitive  data  sets  for  zero  pressure  gradient,  constant 
wall  temperature,  incompressible  turbulent  boundary  layer  flow  over 
smooth  flat  plates  are  those  of  Reynolds,  Kays  and  Kline  1958  [24]. 
In  fact,  these  are  the  only  widely  referenced  data  for  the  condi¬ 
tions  of  interest  and  serve  as  the  basis  for  all  heat  transfer 
correlations  for  these  conditions.  They  are  the  only  data  quoted, 
for  example,  by  Kays  and  Crawford  1980  [31]  and  Rohsenow  and 
Hartnett  1973  [35].  Their  experimental  apparatus  consisted  of  24 


39 


Individually  heated  copper  plates.  The  plate  dimensions  were  6.4 
cm  long  In  the  flow  direction  by  approximately  84  cm  wide.  This 
gave  a  total  surface  length  of  1.5  m.  This  test  surface  was  placed 
in  a  2.3  m  diameter  free-jet  wind  tunnel  with  freestream  turbulence 
intensity  of  between  2  and  5  percent,  depending  on  the  freestream 
velocity.  Stanton  numbers  were  determined  by  measuring  the  power 
input  to  each  plate,  the  plate  temperatures,  and  the  freestream 
recovery  temperature  and  correcting  for  heat  losses. 

Figure  20  shows  a  plot  of  this  data  (with  no  variable  property 
correction)  along  with  the  correlation 

St  -  0. 185 ( 1 og10Rex) ~2* 584 (Pr) ~0*4  (5) 

and  ±5  percent  bands.  A  Prandtl  number  of  0.713  was  used  In  the 
correlation,  which  Is  based  on  the  analogy  StPr0,4  -  C^/2  and  uses 
the  Schultz-Grunow  expression  (equation  (1))  for  Cf/2.  Most  of  the 
data  scatter  within  the  ±5  percent  range.  These  data  represent  8 
individual  runs  with  freestream  velocities  ranging  from  14  m/s  to 
39  m/s  and  with  Rex  up  to  3.5  million. 

Other  smooth  wall  data  sets  chosen  for  comparision  are  from  a 
series  of  experiments  at  Stanford  University.  These  data  sets  are 
from  Moffat  1967  [25],  Kearney  1970  [26]  and  Orlando  1974  [27]. 
These  experiments  were  mainly  concerned  with  the  effects  of  tran¬ 
spiration  on  heat  transfer  In  the  turbulent  boundary  layer.  Their 
surface  was  porous  to  allow  transpiration  and  had  an  rms  roughness 
of  5  micrometers,  which  proved  to  be  aerodynamical ly  smooth  at  the 
low  velocities  used  In  the  experiments.  Each  of  the  workers  took 
baseline  non-transplred  data  to  qualify  the  experimental  apparatus, 
and  it  Is  this  data  which  Is  of  Interest.  All  of  these  experiments 
were  conducted  on  a  test  facility  very  similar  to  the  THTTF.  The 
test  surface  consisted  of  24  Individually  heated  plates.  Each 
plate  was  10.2  cm  long  by  45.7  cm  wide,  resulting  In  a  2.4  m  long 
test  bed.  The  nominal  freestream  turbulence  was  0.7  percent. 
Stanton  numbers  were  determined  by  measuring  the  plate  temperatures 
and  freestream  recovery  temperature  and  the  power  input  to  each 
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plate  and  correcting  for  heat  losses.  Seven  data  sets  have  been 
selected  to  represent  the  early,  mid  and  late  time  periods  of  the 
series.  Figure  21  shows  a  plot  of  the  data  along  with  equation  (5) 
and  the  ±5  percent  bands.  Inspection  of  the  figure  reveals  that 
the  data  all  scatter  within  the  ±5  percent  bands.  These  data 
represent  freestream  velocities  from  7  m/s  to  13.4  m/s  and  with  Rex 
up  to  about  2  million. 

The  comparisons  in  Figures  20  and  21  demonstrate  that  equation 
(5)  is  a  reasonable  representation  of  the  existing  smooth  wall, 
constant  temperature,  zero  pressure  gradient  Stanton  number  data 
and  that  the  data  scatter  within  approximately  ±5  percent  about 
this  correlation.  Therefore,  if  the  THTTF  data  with  their  associ¬ 
ated  uncertainties  of  about  2  to  5%  are  within  the  ±5  percent 
interval  about  equation  (5),  it  can  be  concluded  that  a  successful 
comparison  has  been  achieved  at  the  Nth  order  replication  level 
(Moffat  1988  [36],  Coleman  and  Steele  1989  [37])  and  that  the 
qualification  has  been  proven. 

As  discussed  In  detail  In  Appendix  I,  a  detailed  uncertainty 
analysis  of  the  determination  of  Stanton  numbers  using  the  THTTF 
was  made  during  the  design  and  construction  phase  of  the  test 
facility.  This  analysis  showed  that  the  precision  limits  in  the 
measurements  were  negligible  relative  to  the  bias  limits.  This 
meant  that  replications  of  Stanton  number  at  a  given  experimental 
set  point  (U„)  should  show  negligible  scatter,  since  any  signifi¬ 
cant  errors  were  estimated  to  be  bias  errors  that  would  be  the  same 
in  all  replications. 

This  was  found  to  be  the  case  for  runs  with  U.  -  27,  43,  58 
and  67  m/s  as  Figures  22-25  Illustrate.  The  Stanton  number  results 
for  two  runs  are  plotted  in  each  Instance,  and  it  is  apparent  that 
the  estimate  of  zero  precision  limit  (no  scatter)  is  a  reasonable 
one.  Such  a  conclusion  cannot  be  drawn  for  the  U„  -  12  m/s  condi¬ 
tion,  however.  This  Is  shown  In  Figure  26,  in  which  the  results  of 
eight  separate  Stanton  number  runs  at  U„  -  12  m/s  are  shown.  The 
data  clearly  Indicate  a  run-to-run  scatter  which  cannot  be  termed 
negligible. 
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Investigation  showed  that  at  low  freestream  veloc1t1es--U„  i 
12  m/s  for  the  smooth  wall  tests  and  U„  i  6  m/s  for  the  rough  wall 
tests--the  heat  transfer  coefficients  are  relatively  low  and  the 
time  constant  of  the  THTTF  is  thus  Increased.  At  these  conditions 
the  time  constant  of  the  THTTF  is  large  enough  so  that  the  rela¬ 
tively  long  period  variations  in  facility  line  voltage  to  the  test 
plate  heater  circuits  and  in  the  temperature  of  the  incoming 
make-up  water  for  the  heat  exchanger  loop  affect  the  ability  to 
hold  a  tight  steady  state  condition.  These  annoyances  could  be 
overcome  with  additional  expenditures  for  power  conditioning  equip¬ 
ment  and  a  water  chiller  system;  however,  the  observed  run-to-run 
scatter  in  St  results  is  within  acceptable  limits.  Observations  of 
the  St  results  for  the  nine  runs  shown  in  Figure  26  and  also  three 
U.  «  6  m/s  runs  with  a  rough  wall  produced  a  95%  confidence  esti¬ 
mate  of  a  precision  limit  in  St  of  about  3%  for  these  conditions. 
This  is  present  because  of  system  unsteadiness  and  not  because  of 
measurement  uncertainties.  For  the  U„  -  12  m/s  smooth  wall  data, 
therefore,  there  is  a  3%  precision  limit  contribution  combined  by 
root-sum-square  with  the  estimated  bias  limit  to  obtain  the  overall 
uncertainty  in  St. 

Shown  in  Figure  27  are  the  St  data  from  one  replication  at 
each  freestream  velocity  plotted  along  with  the  95%  confidence 
uncertainty  interval  for  representative  data  points.  These  are 
compared  with  equation  (5)  and  its  ±5%  interval,  and  the  comparison 
shows  the  excellent  agreement  which  is  obtained.  The  data  and  the 
associated  uncertainty  bands  everywhere  overlap  the  correlation 
with  its  associated  ±5%  interval.  This  comparison  validates  the 
qualification  of  the  THTTF  for  Stanton  number  measurements. 

The  present  data  sets  extend  to  a  larger  Reynolds  number  than 

the  previous  data  (Rex  -  1  x  107  versus  3.5  x  106)  and  over  a 

larger  range  of  freestream  velocities.  An  extensive  survey  of  the 

literature  revealed  no  flat  plate  Stanton  number  data  for 
x-Reynolds  numbers  above  3.5  million  for  Incompressible  flows.  In 
the  present  work  the  experimental  Stanton  number  range  has  been 
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essentially  tripled  up  to  x-Reynolds  numbers  of  10  million.  It  can 
be  seen  that  equation  (5)  represents  the  data  well  over  the  entire 
range. 

4.2  TEMPERATURE  PROFILE  COMPARISON 

Mean  temperature  profiles  were  measured  by  traversing  the 
boundary  layer  with  a  thermocouple  probe  as  discussed  In  Appendix 
II.  The  maximum  variation  In  the  freestream  temperature  In  the 
THTTF  from  any  one  temperature  profile  location  to  any  other  tem¬ 
perature  profile  location  was  less  than  0.16°C.  Previous  studies 
have  shown  that  temperature  profiles  along  smooth  flat  plates  of 
constant  temperature  In  a  zero  pressure  gradient,  turbulent  flow 
agree  with  the  law  of  the  wall  for  a  thermal  boundary  layer  as 
given  by  Kays  and  Crawford 

T+  -  2.195  Sin  y+  +  13.2  Pr  -  5.66  (6) 

Figure  28  shows  the  temperature  profiles  at  plate  19  for  U„  -  12 
and  43  m/s  plotted  In  T+  vs.  y+  coordinates  along  with  equation 
(6).  The  data  show  good  agreement  with  the  temperature  "law  of  the 
wall." 
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Stanton  Number  Data  of  Reynolds,  Kays  and  Kline  (1958)  Compared 
the  x-Reynolds  Number  Correlation  (equation  (5)). 


Stanton  Number  Data  of  Moffat  (1967),  Kearney 
Compared  with  the  x-Reynolds  Number  Correlate 


Stanton  Number 
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:igure  27.  Comparison  of  the  THTTF  Smooth  Wall  Stanton  Number  Data  with  the 
x-Reynolds  Number  Correlation  (equation  (5)). 
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SECTION  5 

ADDITIONAL  SMOOTH  WALL  HEAT  TRANSFER  RESULTS 


In  addition  to  the  constant  wall  temperature  qualification 
data  presented  in  the  preceding  sections,  data  were  obtained  for 
several  other  thermal  boundary  conditions.  During  the  search  of 
the  literature  for  qualification  data,  only  one  systematic  experi¬ 
mental  study  of  the  effect  of  thermal  boundary  condition  was  found 
(Reynolds,  Kays,  and  Kline  1958  [24]).  That  work  was  limited  to 
x-Reynolds  numbers  less  than  3,500,000.  Since  the  THTTF  is  able  to 
achieve  x-Reynolds  numbers  of  10,000,000,  the  data  reported  in  this 
section  were  taken  to  add  to  the  existing  data  and  to  extend  the 
x-Reynolds  number  range  to  10,000,000.  In  addition  to  the  constant 
wall  temperature,  five  thermal  boundary  conditions  were  considered: 
(1)  constant  heat  flux,  (2)  step  wall  temperature,  (3)  step  wall 
heat  flux,  (4)  linear  wall  temperature  variation  (ramp  down),  and 
(5)  bi-1 inear  wall  temperature  variation  (ramp  down  then  ramp  up). 

The  two  most  basic  thermal  boundary  conditions  are  those  of 
constant  wall  temperature  and  constant  heat  flux.  These  two  cases 
are  very  good  approximations  for  many  real  heat  transfer  problems, 
are  relatively  easy  to  construct  in  the  laboratory,  and  allow  ready 
solution  of  the  boundary  layer  equations.  Therefore,  they  are 
often  used  in  test  cases  for  theoretical  heat  transfer  computa¬ 
tions. 

The  step  wall  temperature  case  is  one  of  the  fundamental 
problems  of  convective  heat  transfer.  Under  the  assumptions  of 
incompressible  flow  with  constant  fluid  properties,  the  momentum 
and  energy  equations  become  uncoupled,  and  the  energy  equation 
becomes  linear.  Therefore,  the  problem  of  heat  transfer  In  the 
boundary  layer  with  arbitrary  thermal  boundary  condition  becomes 
amenable  to  solution  by  superposition.  The  simplest  boundary 
condition  for  which  the  solution  can  serve  as  a  kernel  function  for 
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superposition  Is  the  step  wall  temperature  case.  Thus,  super¬ 
position  Integral  methods  can  be  formulated  based  on  the  step  wall 
temperature  results  to  predict  the  heat  transfer  (Stanton  number) 
for  arbitrary  wall  temperature  or  heat  flux. 

In  the  following,  the  theory  of  heat  transfer  In  the  turbulent 
incompressible  smooth  wall  boundary  layer  Is  briefly  reviewed  and 
the  experimental  results  are  presented,  discussed,  and  compared. 

5.1  THEORY 

The  theoretical  treatment  of  heat  transfer  in  turbulent  incom¬ 
pressible  smooth  flat  plate  boundary  layer  flow  is  mature  and  well 
documented  (Cebeci  and  Bradshaw  1984  [38],  and  (Cays  and  Crawford 
1980  [31]).  Here  the  theory  is  divided  into  3  subtopics:  analo¬ 
gies,  solutions  to  the  integral  boundary  layer  equations,  and 
numerical  solutions  of  the  time  averaged  differential  boundary 
layer  equations. 

5.1.1  Analogies 

Reynolds  analogy  (Reynolds  1894  [39])  was  one  of  the  first 
theoretical  treatments  of  turbulent  heat  transfer.  For  the  case  of 
constant  freestream  velocity  and  temperature,  constant  wall  tem¬ 
perature,  and  Pr  -  1,  the  similarity  between  the  x-momentum  equa¬ 
tion  and  the  energy  equation  yields 

St  *  Cf/2  (7) 

The  restrictions  on  equation  (7)  prohibit  its  wide  application. 
Particularly  troublesome  is  the  requirement  that  Pr  -  1.  The  two 
most  famous  extensions  of  the  analogy  for  Pr  of  order  1  or  greater 
are  the  one  of  Von  Karman  1939  [40]  and  the  one  of  Colburn  1933 
[41]. 

Von  Karman  derived  his  analogy  by  integrating  the  thermal  and 
velocity  laws-of-the-wall  in  a  three  layer  model  to  yield 
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(8) 


’  §7 +  5  ^1/2  {Pr ' 1  +  lnC1  +  5(Pr  ‘ 1)/6]1 

Kays  and  Crawford  1980  [31]  have  shown  that  equation  (8)  can  be 
approximated  very  well  for  gases  (0.5  <  Pr  <  1)  by 

St  Pr0*4  -  Cf/2  (9) 

Colburn  based  his  analogy  completely  on  dimensional  analysis 
and  empirical  considerations.  His  analogy  is 

St  Pr0*66  -  Cf/2  (10) 

Although  equation  (10)  was  first  determined  by  Colburn  purely  from 
empirical  considerations,  it  has  been  deduced  from  turbulent  scal¬ 
ing  laws  for  Pr  >  1  (Bejan  1984  [42]  and  Arpaci  and  Larsen  1984 
[43]). 

Other  analogies  can  be  derived  from  the  velocity  and  tempera¬ 
ture  laws-of-the-wall  (Kader  and  Yaglom  1972  [44],  for  example). 

All  of  these  analogies  are  still  limited  to  the  constant  wall  tem¬ 
perature  boundary  condition.  However,  this  limitation  is  often  not 
stressed  and  the  analogies  are  often  compared  with  constant  heat 
flux  experiments. 

5.1.2  Integral  Equations 

For  the  Integral  analysis,  we  assume  constant  freestream 
velocity  and  temperature  and  a  constant  property  fluid  flow  over  a 
smooth  flat  plate.  For  arbitrary  wall  temperature  or  heat  flux,  we 
follow  the  procedure  of  Reynolds,  Kays  and  Kline  1958  [24]  as  pre¬ 
sented  by  Kays  and  Crawford  1980  [31].  The  procedure  is  to  use  the 
solution  of  the  Integral  boundary  layer  equations  for  the  step  wall 
temperature  with  an  unheated  starting  length  as  the  kernel  function 
in  a  superposition  integral.  Using  the  1/7  power  law  approximation 
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of  the  velocity  and  temperature  profiles,  Reynolds,  Kays,  and  Kline 
established  that  for  an  unheated  starting  length,  £,  the  local  heat 
transfer  coefficient  can  be  expressed  as 


St(£;x)  -  -  Stt(x)[l  -  (C/x)9/10]“1/9  (11) 

PU.Cp 

where  Stt (x)  is  the  constant  wall  temperature  Stanton  number  given 
by  the  analogy  of  choice  from  equations  (7)  -  (10).  The  heat  flux 
for  a  variable  wall  temperature,  Tw(x),  is 

rx  d(T  -  T  1  N 

q"w(x)  -  JQh(C;x)  -  —  d$  +  Ih(C1;x)A(Tw,i  -  TJ  (12) 

where  A(TW j  -  T„)  is  the  1 ‘ th  finite  step  In  wall  temperature,  N 
is  the  number  of  steps,  and  the  ^  are  the  locations  of  the  steps. 

For  a  piecewise  linear  wall  temperature  with  a  finite  number 
of  steps 

M  M 

Tw  -  T.  -  l  mi(x  -  a,)  +  l  b1  (13) 

j-1  j-1 

application  of  equation  (12)  with  equation  (11)  yields 

St(x)  J  mj  P  (8/9,10/9) 

9{'w  -  T„)  j,!  J  j 

♦  7s— |  Z  b,ci  -  (e,/x)9/10r1/9  (i4) 

where  0r  (8/9,10/9)  Is  the  Incomplete  beta  function  with 
rj  -  1  -  (aj/x)9/10 

pr(a,b)  -  za-1(l-z)b"1  dz  (15) 
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For  the  step  heat  flux  case  with  unheated  length  t, 

St  T(l/9)  T(8/9) 

- _J - ' — i -  (16) 

$tt  8r(l/9,  10/9) 

with  r  -  1  -  (<fr/x)°*9 

5.1.3  Differential  Equations 

The  finite  difference  solution  of  the  Incompressible  boundary 
layer  equations  is  now  routine  for  smooth  wall  cases.  Cebeci  and 
Bradshaw  1984  [38]  present  a  complete  discussion  of  these  solu¬ 
tions.  The  solutions  presented  in  this  section  are  based  on  a 
mixing  length  turbulence  model  with  van  Driest  damping  and  a  con¬ 
stant  turbulent  Prandtl  number,  Prt  -  0.9.  For  these  computations, 
the  BLACOMP  code  as  verified  by  Gatlin  1983  [45]  was  used.  The 
particular  details  of  the  solutions  presented  here  can  be  found  in 
Love  et  al.  1988  [46].  As  was  the  case  with  the  experimental  data 
reduction,  all  fluid  properties  were  assumed  to  be  constant  and 
were  evaluated  at  the  freestream  static  temperature. 

5.2  EXPERIMENTAL  RESULTS 

Stanton  number  data  were  taken  for  both  constant  wall  tempera¬ 
ture  and  constant  heat  flux  boundary  conditions  for  nominal  free¬ 
stream  velocities  of  12,  27,  42,  and  67  m/s.  The  experiments  were 
conducted  so  that  the  constant  heat  flux  cases  had  approximately 
the  same  wall  temperature  far  from  the  boundary  layer  origin  as  the 
corresponding  constant  wall  temperature  cases.  These  data  are 
compared  with  the  analogies,  the  integral  solutions,  and  the  nu¬ 
merical  solutions  of  the  boundary  layer  equations.  However,  the 
primary  comparison  Is  the  comparison  of  the  data  sets  with  each 
other. 

Figure  29  presents  a  plot  of  St  versus  Rex  for  the  two  bound¬ 
ary  conditions  (constant  q"w  and  Tw).  The  solid  symbols  are  the 
results  for  the  constant  heat  flux  boundary  condition  and  the  open 
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symbols  are  for  the  constant  wall  temperature  boundary  condition. 
The  solid  curve  Is  the  analogy  of  equation  (9),  StPr0,4  -  Cf/2, 
where  the  skin  friction  coefficient  Is  determined  using  the 
Schultz-Grunow  1941  [  29]  correlation 

Cf/2  -  0.185[Log10(Rex)  ]'2*584  (17) 

The  figure  shows  that  the  analogy  fits  the  constant  wall  tempera¬ 
ture  data  very  well  for  the  Reynolds  number  range  (100,000  to 
10,000,000).  At  the  maximum  deviation  the  correlation  is  about  5 
percent  too  high.  The  dashed  curve  Is  based  on  the  Integral  solu¬ 
tion,  equation  (16),  for  the  case  of  constant  q“w 


r(i/9)  r(8/9)  :t 
P1(l/9,10/9)  t 


1.043  Stt 


(18) 


Far  from  the  boundary  layer  origin,  equation  (18)  agrees  almost 
exactly  with  the  data.  For  Reynolds  numbers  between  3,000,000  and 
10,000,000,  the  constant  wall  heat  flux  data  are  consistently  4  to 
5  percent  greater  than  the  constant  wall  temperature  data.  How¬ 
ever,  near  the  origin  of  the  boundary  layer  the  constant  wall  heat 
flux  data  are  10  to  15  percent  above  the  constant  wall  temperature 
data. 

The  two  boundary  conditions  are  compared  more  directly  in 
Figure  30.  The  ratio  of  constant  wall  heat  flux  to  constant  wall 
temperature  Stanton  numbers  Is  plotted  directly.  The  dashed  curve 
represents  the  Integral  solution,  equation  (18),  and  the  solid 
lines  represent  the  numerical  solutions  of  the  boundary  layer 
equations.  For  Reynolds  numbers  greater  than  3,000,000,  the  data 
agree  very  closely  with  equation  (18).  But,  for  the  low  Reynolds 
numbers,  the  data  Indicate  a  ratio  of  1.10  to  1.15  instead  of  1.04. 
The  numerical  solutions  agree  with  the  data  for  all  of  the  Reynolds 
numbers.  The  numerical  solutions  Indicate  a  slight  unit  Reynolds 
number  effect  which  Is  perhaps  suggested  by  the  data  but  not 
proven. 
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As  previously  discussed,  the  Reynolds  analogies  strictly  apply 
only  for  the  constant  wall  temperature  boundary  condition.  How¬ 
ever,  they  are  also  considered  to  apply  approximately  for  the 
constant  heat  flux  boundary  condition.  Constant  wall  heat  flux 
data  are  often  reported  in  terms  of  the  Reynolds  analogy  factor, 
2St/Cf,  (for  example,  Subramanlan  and  Antonia  1981  [47],  and 
Simonich  and  Bradshaw  1978  [48]).  In  Figure  31  the  data  are  pre¬ 
sented  in  analogy  coordinates,  2St/Cf,  directly.  Again  the  solid 
symbols  are  the  constant  q"w  results  and  the  open  symbols  are  the 
constant  Tw  results.  The  skin  friction  coefficients  are  the  values 
measured  with  a  Preston  tube.  The  dashed  lines  are,  as  indicated, 
the  Colburn  analogy,  Pr-0*^6,  and  the  approximation  of  Von  ICarman' s 
analogy,  Pr-0*4.  The  solid  lines  represent  the  results  from  the 
numerical  solutions.  The  analogy  factors  for  the  constant  Tw  data 
are  more  or  less  constant  for  the  whole  Reynolds  number  range.  The 
data  scatter  about  the  value  of  about  1.2  which  is  approximately 
halfway  between  the  two  analogies.  The  numerical  solutions  for  the 
constant  Tw  case  are  in  very  good  agreement  with  the  data.  They 
are  practically  constant  at  2St/Cf  -  l.i8  to  1.2  and  show  only  a 
small  unit  Reynolds  number  effect.  The  analogy  factors  for  the 
constant  q"w  boundary  condition  show  a  Reynolds  number  dependence. 
The  ratio  decreases  from  about  1.4  at  low  Reynolds  numbers  to  about 
1.24  at  the  higher  Reynolds  numbers.  The  numerical  solutions  are 
in  good  agreement  with  the  data  and  show  a  small  unit  Reynolds 
number  effect.  Simonich  and  Bradshaw  1978  [48]  and  Subramanlan  and 
Antonia  1981  [47]  presented  limited  heat  transfer  data  for  turbu¬ 
lent  boundary  layers  with  constant  q"w  boundary  conditions. 

Simonich  and  Bradshaw  found  that  2St/Cf  varied  between  1.3  at  low 
Reynolds  numbers  to  1.2  at  Rex  -  3,500,000.  Subramanlan  and 
Antonia  found  2St/Cf  to  vary  from  1.5  at  low  Rex  to  1.4  at  Rex  - 
3,000,000.  These  results  are  In  substantial  agreement  with  the 
present  results. 

For  unheated  starting  length  boundary  conditions,  Stanton 
number  measurements  were  made  at  nominal  freestream  velocities  of 
28  m/s  and  67  m/s.  Three  cases  were  run  at  each  velocity  for  each 
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boundary  condition  (step  Tw  and  q"w).  The  lengths  of  the  unheated 
regions  were  chosen  so  that  an  appropriate  spread  In  Reynolds 
numbers,  Re^,  was  obtained.  At  U„  -  28  m/s,  the  unheated  starting 
lengths  were  0.3  m,  0.7  m,  and  1.3  m.  At  U.  -  67  m/s,  they  were 
0.5  m,  0.8  m,  and  1.3  m. 

Figure  32  shows  a  summary  of  the  Stanton  number  data  for  a 
constant  wall  temperature  boundary  condition  and  the  step  Tw  cases 
for  U_  -  28  m/s  and  67  m/s.  The  first  heated  plate  Is  highlighted 
in  each  case  by  plotting  its  data  as  a  solid  symbol.  Data  from  the 
last  plate  is  not  plotted  for  any  case.  The  figure  shows  that  a 
step  in  wall  temperature  has  a  large  effect  on  the  Stanton  number 
in  the  heated  region  near  the  step.  But  as  the  thermal  boundary 
layer  develops,  the  Stanton  numbers  approach  the  results  for  the 
constant  wall  temperature  boundary  layers.  The  starting  lengths 
were  chosen  so  that  the  last  case  at  U„  ■  28  m/s  and  the  first  case 
at  U.  -  67  m/s  had  approximately  the  same  value  of  Re^.  Based  on 
the  data  of  Reynolds  et  al.,  the  results  of  these  two  cases  should 
coincide  when  plotted  as  St  versus  Rex.  The  figure  shows  that  this 
is  true  for  the  present  results. 

Figure  33  shows  a  comparison  of  the  results  of  the  step  Tw 
experiments  with  the  integral  solution  in  equation  (11),  dashed 
lines,  and  with  the  finite  difference  solutions,  solid  lines,  for 
U„  -  28  m/s.  Figure  34  shows  the  same  comparison  for  U_  -  67  m/s. 
The  results  are  presented  in  terms  of  St/Stt  for  a  direct  compari¬ 
son  with  equation  (11).  The  Stc  data  were  used  to  normalize  the  St 
data,  and  the  finite  difference  solutions  for  constant  wall  tem¬ 
perature,  Stt,  were  used  to  normalize  the  finite  difference  solu¬ 
tions,  St.  The  figures  show  that  both  the  finite  difference  and 
Integral  solutions  are  In  good  agreement  with  the  data  in  all 
cases.  The  Integral  solutions  are  consistently  low  by  a  small 
amount  In  the  region  of  the  step  In  wall  temperature.  From  the 
comparisons  In  the  figures.  It  can  be  concluded  that  equation  (11) 
is  still  appropriate  for  values  of  Rex  -*  107. 
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Figure  35  shows  a  summary  of  the  Stanton  number  data  for  a 
constant  heat  flux  boundary  condition  and  the  step  q"w  cases  for  U. 

-  28  m/s  and  67  m/s.  The  figure  shows  that  as  the  thermal  boundary 
layer  develops,  the  unheated  starting  length  Stanton  numbers  ap¬ 
proach  the  results  for  the  constant  heat  flux  boundary  condition. 

The  first  heated  plate  is  highlighted  in  each  case  by  plotting  its 
data  as  a  solid  symbol.  Data  from  the  last  plate  is  not  plotted 
for  any  case. 

Figure  36  shows  a  comparison  of  the  results  of  the  experiments 
with  the  integral  solution  in  equation  (16),  dashed  lines,  and  with 
the  finite  difference  solutions,  solid  lines,  for  U_  -  28  m/s. 

Figure  37  shows  the  same  comparison  for  U.  «  67  m/s.  As  before, 
the  results  are  presented  in  terms  of  St/St^  for  a  direct  compari¬ 
son  with  equation  (16).  The  figures  show  that  the  finite  difference 
solutions  are  in  very  good  agreement  with  the  data  In  all  cases. 

The  Integral  solutions  are  also  in  reasonable  agreement  with  the 
data,  with  the  maximum  difference  between  the  data  and  the  Integral 
solutions  being  about  10  percent. 

As  examples  of  variable  wall  temperature  boundary  conditions, 
one  linearly  decreasing  (ramp  down)  and  two  bi-1 Inear  (ramp  down 
then  ramp  up)  wall  temperature  distributions  were  used  for  experi¬ 
ments  at  freestream  velocities  of  28  and  67  m/s.  The  measured 
axial  wall  temperature  distributions  are  shown  in  Figure  38.  A 
discussion  of  the  complete  data  set  is  given  in  Love  et  al.  1988 
[46].  Only  two  sample  cases  are  given  here. 

Figure  39  displays  the  data  for  the  linear  wall  temperature 
distribution  for  a  freestream  velocity  of  67  m/s.  The  dashed  line 
represents  the  Integral  solutions  of  equation  (14),  and  the  solid 
line  represents  the  numerical  solutions.  The  solid  symbols  Indi¬ 
cate  the  locations  of  changes  In  the  slope  of  the  wall  temperature 
distribution.  Figure  40  shows  the  same  information  for  the 
bl -1  Inear  wall  temperature  distribution  with  1.0°C  temperature 
difference  from  plate  to  plate.  Similar  to  the  unheated  starting 
length  cases,  the  ratio  $t/Stt  is  used  for  a  direct  comparison  with 
the  results  of  the  superposition  Integral.  Although  the  numerical 
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solutions  predict  a  greater  response  to  changes  In  the  wall  tem¬ 
perature,  both  the  numerical  and  Integral  solutions  are  In  reason¬ 
able  agreement  with  the  experimental  data. 

5.3  SUMMARY 

This  study  extended  the  x-Reynolds  number  range  of  such  data 
for  incompressible  smooth  flat  plate  flows  from  a  maximum  of 
3,500,000  to  10,000,000.  As  seen  in  Section  4,  the  Stanton  numbers 
for  constant  were  in  good  agreement  with  the  Von  Karman  analogy 
over  the  entire  Reynolds  number  range  when  the  skin  friction  coef¬ 
ficient  given  by  the  Schultz-Grunow  correlation  is  used.  Stanton 
numbers  for  the  constant  heat  flux  boundary  condition  were  5  to  15% 
higher  than  those  for  the  constant  wall  temperature  boundary  condi¬ 
tion.  The  constant  wall  temperature  experiments  yielded  Reynolds 
analogy  factor,  2St/Cf,  data  which  scattered  about  the  constant 
1.2.  The  analogy  factor  decreased  from  about  1.4  at  lower  Reynolds 
numbers  to  1.2  downstream  for  the  constant  heat  flux  boundary 
condition.  The  numerical  solutions  were  in  good  agreement  with  the 
data  for  both  the  constant  wall  temperature  and  constant  heat  flux 
cases. 

The  integral  solutions  and  numerical  solutions  were  in  sub¬ 
stantial  agreement  with  the  experimental  Stanton  number  data  for 
the  constant  wall  temperature  and  constant  heat  flux  boundary 
conditions  with  unheated  starting  lengths.  The  integral  solutions 
and  the  numerical  solutions  gave  good  agreement  with  the  experimen¬ 
tal  Stanton  number  data  for  the  linear  and  bi-1 inear  wall  tempera¬ 
ture  distributions.  The  effect  of  a  decreasing  wall  temperature  in 
the  flow  direction  was  to  lower  the  local  Stanton  number  to  less 
than  the  constant  wall  temperature  Stanton  number  Stt.  As  the  wall 
temperature  Increased  in  the  flow  direction,  the  local  Stanton 
number  became  greater  than  Stt. 
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Figure  29.  Stanton  Number  Results  for  Constant  Wall  Temperature  (open  symbols) 
and  Constant  Wall  Heat  Flux  (solid  symbols). 
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gure  30.  Direct  Comparison  of  Constant  Heat  Flux  and  Constant  Wall  Temperature  Thermal 
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Figure  33.  Comparison  of  the  Results  of  the  Step  T  Experiments  with  the  Integral  and 
Finite  Difference  Solutions.  w 
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Stanton  Numbers  for  Constant  q  "  with  Unheated  Starting  Lengths 
Filled  Symbols  Indicate  First  wHeated  Plate. 
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Figure  37.  Comparison  of  Experimental  Stanton  Number  Ratios  for  Constant  q  "  with  Unheated 
Starting  Lengths  with  Integral  Solutions  and  Numerical  Solution£--U  =  67  m/s. 
Filled  Symbols  Indicate  First  Heated  Plate. 
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Figure  38.  Experimental  Linear  and  Bi-linear  Wall  Temperature  Distributions. 
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Figure  39.  Comparison  of  Experimental  Stanton  Number  Ratio  for  the  Linear  Wall  Temperature 
Distribution  with  Integral  and  Numerical  Solutions--!^  =  67  m/s.  Filled  Symbol 
Indicates  Change  in  Slope  of  T  . 
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Figure  40.  Comparison  of  Experimental  Stanton  Number  Ratio  for  Bi-1 inear  Wall  Temperature 
Distribution  with  1°C  Plate  to  Plate  Temperature  Difference  with  Integral  and 
Numberical  Solutions--!!^  =  67  m/s.  Filled  Symbol  Indicates  Change  in  Slope  of  T 


SECTION  6 

ROUGH  WALL  RESULTS 


Stanton  number  data  have  been  obtained  on  one  rough  surface 
for  freestream  velocities  of  6,  12,  28,  43,  58  and  67  m/s.  This 
surface  is  the  one  roughened  with  1.27  mm  diameter  hemispherical 
elements  spaced  4  diameters  apart  in  a  staggered  array  as  shown  in 
Figure  41.  In  addition,  skin  friction  coefficient  distributions 
have  been  obtained  at  U.  -  12  and  58  m/s.  Profiles  of  mean  veloc¬ 
ity,  mean  temperature,  and  the  Reynolds  stress  quantities  u'2,  v'2, 
w'2,  and  u'v1  have  also  been  measured.  These  data  are  presented 
and  discussed  in  this  section  along  with  a  comparison  with  Stanton 
number  and  skin  friction  coefficient  predictions  made  using  the 
discrete  element  roughness  model  of  Taylor,  Coleman  and  Hodge  1984 
[49]. 

6.1  REVIEW  OF  THE  DISCRETE  ELEMENT  PREDICTION  MODEL 

The  development  and  calibration  of  the  discrete  element  pre¬ 
diction  model  has  previously  been  described.  Only  an  overview  is 
given  here.  The  basic  idea  of  the  discrete  element  approach  is  to 
treat  the  roughness  as  a  collection  of  individual  entities  and  to 
account  for  the  blockage,  form  drag,  and  heat  transfer  on  the 
elements.  The  discrete  element  model  presented  here  is  formulated 
for  roughness  elements  with  three-dimensional  shapes  (as  opposed  to 
transverse  ribs,  for  example)  for  which  the  element  cross-section 
can  be  approximated  as  circular  at  every  height,  y. 

The  differential  equations  Including  roughness  effects  are 
derived  by  applying  the  basic  conservation  statements  for  mass, 
momentum  and  energy  to  a  control  volume  (CV)  such  as  that  shown  in 
Figure  42.  Upon  Incorporating  the  boundary  layer  assumptions,  the 
continuity,  momentum  and  energy  equations  for  a  steady 
(Reynolds-averaged),  two-dimensional  turbulent  boundary  layer  are 


75 
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and 
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Examination  of  equations  (20)  and  (21)  shows  that  empirical  models 
for  -pu 1 v 1 ,  -pv'h',  and  the  roughness  element  drag  coefficient  CD 
and  Nusselt  number  Nud  are  necessary  for  closure. 

The  blockage  parameters  px  and  py  and  the  element  shape 
descriptor  d(y)  require  no  empirical  fluid  mechanics  Input  as  they 
are  determined  solely  from  the  c’.ometry  of  the  rough  surface. 
Taylor  et  al.  1984  [49]  have  shown  for  uniform  three-dimensional 
roughness  elements  with  circular  cross-section  that 
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The  boundary  conditions  for  the  discrete  element  approach  for 
rough  wall  flows  are  Identical  to  those  for  smooth  wall  flows.  The 
wall  location  (y  -  0)  Is  the  smooth  surface  on  which  the  roughness 
elements  occur.  At  y  «  0,  u-v-0  and  H  -  Hw.  As  y  -*  -,  u  ->  U. 
and  H  -  H„. 

The  "wall  shear  stress"  Is  defined  as  the  sum  of  the  shear  and 
drag  forces  on  the  wall  In  the  mean  flow  direction  divided  by  the 
plan  area  of  the  wall.  The  corresponding  skin  friction  coefficient 
is  then 
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and  the  Stanton  Number  is 
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These  definitions  for  and  St  can  be  formulated  from  physical 
reasoning.  However,  they  also  arise  naturally  from  equations  (19) 
-(21)  In  the  formulation  of  the  integral  boundary  layer  equations 
using  the  discrete  element  model. 

In  order  to  solve  equations  (19),  (20)  and  (21),  turbulence 
models  for  -pu'v'  and  -pv'h1  and  roughness  models  for  CD  and  Nud 
are  required.  Because  of  Its  wide  acceptance  and  proven  predictive 
capability  for  boundary  layer  flows  over  smooth  surfaces,  the 
Prandtl  mixing  length  model  with  van  Driest  damping  and  a  constant 
turbulent  Prandtl  number  is  used  for  turbulence  closure.  Thus 
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where 


-  0.40y[l  -  exp(-y+/26)]  ;  <  0.096  (26) 

lm  -  0.096;  otherwise,  (27) 

and 

Prt  -  0.9  (28) 

The  authors,  as  did  Lin  and  Bywater  1980  [50],  chose  to  formu¬ 
late  the  CD  and  Nud  models  as  functions  of  the  local  element 
Reynolds  number 

Red  .  u(y)d(y)/v  (29) 

which  includes  roughness  element  size  and  shape  information  through 
d(y).  The  general  shape  of  the  drag  coefficient  and  Nusselt  number 
versus  Reynolds  number  curves  of  Zukauskas  1972  [51]  were  used  as 
starting  points  to  determine  the  final  expressions  for  CD  and  Nud 
through  comparison  with  calibration  data  sets  from  well-defined 
rough  surfaces.  As  discussed  in  Taylor  et  al.  1984  [49],  the  model 
which  gave  the  best  overall  agreement  was 


log  CD  -  -0.125  log  (Red) 
+  0.375 

Cq  ■  0.6 

Nud  -  2.475  Red°*4Pr0*36; 
Nud  -  1.043  Red0,5Pr0,37; 
Nud  -  0.963  Red°-6Pr0*36; 


Red  <  6  x  104 

(30) 

Red  l  6  x  104 
Red  <  100 

100  <  Red  <  1000  (31) 

1000  <  Red  <  200,000 
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The  C0  model  has  been  tested  for  values  of  Red  up  to  about  25,000 
(Taylor,  et  al.  1984  [49])  and  up  to  about  17,000  (Scaggs  et  al. 
1988  [21]).  The  Nud  model  has  been  tested  for  values  of  Red  up  to 
about  1000  (Taylor  et  al.  1984  [49])  and  has  been  extended  up  to 
about  2400  in  this  report. 

6.2  HEAT  TRANSFER  RESULTS 

One  of  the  most  important  aspects  of  this  research  program  is 
the  development  of  boundary  layer  heat  transfer  data  which  can  be 
used  to  formulate  and  refine  predictive  models.  Stanton  number 
results  have  been  obtained  for  freestream  velocities  of  6,  12,  28, 
43,  58,  and  67  m/s  on  the  rough  surface  composed  of  1.27  mm  diame¬ 
ter  hemispheres  spaced  in  a  staggered  array  4  diameters  apart.  In 
addition,  mean  temperature  profiles  have  been  obtained  at  U.  »  12 
m/s. 

Figure  43  shows  a  composite  plot  of  the  Stanton  number  data 
for  this  surface.  The  error  bars  on  selected  points  indicate  the 
uncertainties  as  computed  using  the  techniques  discussed  in  Appen¬ 
dix  I.  The  curve  represents  smooth  wall  results  as  given  by  equa¬ 
tion  (5).  Figure  44  shows  the  Stanton  number  data  plotted  versus 
A2/k  where  A2  is  the  enthalpy  thickness  and  k  is  the  roughness 
height.  Pimenta  1975  [13]  found  that  for  the  Stanford  surface  (1.27 
mm  diameter  spheres  packed  in  the  most  dense  array),  Stanton  number 
data  plotted  in  these  coordinates  fell  on  a  single  curve  for  bound¬ 
ary  layers  In  the  fully  rough  regime.  Inspection  of  the  figure 
shows  that  as  the  freestream  velocity  Increases  the  data  approach  a 
single  curve,  with  the  data  for  U„  ■  58  m/s  and  67  m/s  essentially 
coinciding.  Based  on  this  criteria  the  boundary  layer  on  this 
surface  would  be  considered  to  be  In  the  transitionally  rough 
regime  for  the  lower  freestream  velocities  and  in  the  fully  rough 
regime  for  U„  -  58  m/s  and  67  m/s.  More  discussion  of  the  transi¬ 
tionally  rough  and  fully  rough  regimes  will  be  given  later  when  the 
measurements  in  the  momentum  boundary  layer  are  discussed. 
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In  Figures  45-50,  the  Stanton  number  data  sets  are  compared 
with  predictions  made  with  the  discrete  element  method  discussed 
previously.  That  model  was  used  without  any  modification;  there¬ 
fore,  these  computations  are  true  predictions.  For  the  lowest 
velocity,  U.  -  6  m/s,  the  computations  underpredict  the  data.  For 
the  intermediate  velocities— 12,  28,  and  43  m/s— the  predictions 
agree  almost  exactly  with  the  data.  For  the  higher  velocities,  U. 

*  58  and  67  m/s,  the  computations  overpredict  the  data,  but  only  by 
about  5  to  8%. 

Figure  51  shows  a  plot  of  a  typical  temperature  profile  for  U.. 

*  12  m/s  in  T+  versus  y+  coordinates.  The  symbols  represent  the 
measurements  over  the  rough  surface  and  the  solid  line  represents 
the  temperature  "law-of-the-wal  1  “  for  smooth  surfaces 

T+  -  2.195  JLn  y+  +  13.2  Pr  -  5.66  (32) 

The  comparison  indicates,  as  expected,  a  downward  shift  In  the  tem¬ 
perature  profile  for  the  rough  wall  data. 

6.3  FLUID  DYNAMICS  RESULTS 

Fluid  dynamics  data  have  been  obtained  in  the  boundary  layer 
for  freestream  velocities  of  12  and  58  m/s.  The  data  were  obtained 
with  the  hot  wire  probes  using  the  techniques  discussed  in 
Appendices  II  and  III.  The  data  include  distributions  of  skin 
friction  coefficient  along  the  surface  and  profiles  of  mean  veloc¬ 
ity  and  Reynolds  stress  quantities. 

Figure  52  shows  a  plot  of  skin  friction  coefficient  versus 
Reynolds  number  based  on  momentum  thickness,  Re^,  for  ■  12  and 
58  m/s.  The  error  bars  represent  the  estimated  ±  10%  uncertainty 
in  Cf.  The  curves  represent  the  predictions  using  the  discrete 
element  model  and  the  smooth  wall  correlation,  equation  (3). 
Comparison  of  the  dat?  and  predictions  shows  that  the  agreement  is 
excellent,  with  the  predictions  agreeing  with  the  data  within  the 
uncertainty.  Such  agreement  is  not  surprising,  since  the  roughness 
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element  drag  coefficient  model  In  the  discrete  element  approach  has 
been  verified  using  a  very  wide  range  of  data  (Scaggs,  Taylor  and 
Coleman  1988  [21]). 

Figure  53  shows  a  typical  mean  velocity  profile  in  u+  versus 
y+  coordinates  for  U„  -  12  m/s.  The  symbols  represent  the  data  and 
the  line  represents  the  usual  smooth  wall  “law-of-the-wal 1 " 

u+  -  2.5  In  y+  +  5.5  (33) 

As  In  the  case  of  the  temperature  profile,  the  expected  downward 
shift  from  the  smooth  wall  case  is  observed. 

In  addition  to  the  mean  velocity  profiles,  measurements  of 
profiles  of  the  Reynolds  stress  quantities  u1^,  v'^,  w'^,  and  u'v1 
were  made.  Figure  54  shows  the  profile  of  (u'^J^/u*  versus  y/S 
for  U.  -  12  m/s.  The  profile  shows  a  behavior  which  is  very  simi¬ 
lar  to  the  profile  over  a  smooth  surface.  Figure  55  shows  a  com¬ 
posite  plot  for  all  three  Reynolds  normal  stress  components.  The 
behavior  observed  in  these  two  figures  indicates  that  the  peak 
value  of  turbulent  kinetic  energy  occurs  very  close  to  the  wall. 
This  is  a  strong  indication  that  the  boundary  layer  is  in  the 
transitionally  rough  regime,  Pimenta  1975  [13].  Figure  56  shows 
the  measured  Reynolds  shear  stress  profile,  which  behaves  as  ex¬ 
pected.  Figure  57  shows  a  plot  of  (u'^)^/u*  versus  y/$  for  = 
58  m/s.  The  figure  shows  a  profile  that  is  typical  of  fully  rough 
boundary  layers,  Pimenta  1975  [13].  The  peak  in  the  turbulence 
intensity  has  moved  to  a  distance  that  is  further  from  the  wall  and 
is  much  broader. 

A  measure  of  the  state  of  the  boundary  layer  flow  (aero¬ 
dynamical  ly  smooth,  transitionally  rough,  or  fully  rough)  which  was 
proposed  by  Taylor  et  al.  1984  [49]  is  the  ratio 

Rt  -  —  (34) 

xtot 
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where  xr  Is  the  part  of  the  apparent  total  surface  shear  stress, 
xtot,  caused  by  the  drag  on  the  roughness  elements  as  computed 
using  the  discrete  element  method.  Scaggs  et  al.  1988  [21]  sug¬ 
gested  that  values  of  Rx  between  about  0.1  and  0.6  Indicate  a 
transitionally  rough  flow.  For  a  freestream  velocity  of  12  m/s  the 
computed  values  of  Rx  are  about  0.4,  Indicating  a  transitionally 
rough  flow.  For  U„  -  58  and  67  m/s  Rx  -  0.6,  indicating  that  these 
flows  should  be  approaching  fully  rough  behavior.  These  conclu¬ 
sions  are  consistent  with  the  observations  of  the  Stanton  number 
data  and  the  turbulence  profiles. 
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Figure  41.  Surface  Roughness  Geometry. 
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Figure  44.  Composite  Plot  of  Rough  Wall  Stanton  Number  Data  vs.  Enthalpy  Thickness 
Normalized  by  the  Roughness  Element  Height. 


Rough  Wall  Stanton  Number  Data  at  a  Freestream  Velocity  of  6  m/s  Compared 
with  the  Prediction  from  the  Discrete  Element  Method. 


Stanton  Number  Data  at  a  Freestream  Velocity  of  12  m/s  Compared 
ediction  from  the  Discrete  Element  Method. 


Stanton  Number  Data  at  a  Freestream  Velocity  of  28  m/s  Compared 
ediction  from  the  Discrete  Element  Method. 


(S/Ul) 


Stanton  Number  Data  at  a  Freestream  Velocity  of  43  m/s  Compared 
ediction  from  the  Discrete  Element  Method. 
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49.  Rough  Wall  Stanton  Number  Data  at  a  Freestream  Velocity  of  58  m/s  Compared 
with  the  Prediction  from  the  Discrete  Element  Method. 


Discrete  Element  Method. 
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Figure  54.  Rough  Wall  Axial  Turbulence  Intensity  Profile  for  a  Freestream  Velocity  of 
12  m/s. 
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Figure  55.  Rough  Wall  Profiles  of  the  Three  Normal  Reynolds  Stresses. 


Figure  56.  Rough  Wall  Reynolds  Shear  Stress  Profile. 


SECTION  7 
SUMMARY 


In  this  report,  the  uncertainty  analysis,  design,  and 
qualification  of  the  Turbulent  Heat  Transfer  Test  Facility  (THTTF) 
with  a  smooth  wall  test  section  have  been  presented,  and  the 
initial  results  for  a  wall  roughened  with  hemispherical  elements 
spaced  four  base  diameters  apart  have  been  discussed  and  compared 
with  calculations  from  the  discrete  element  prediction  method. 

A  detailed  description  of  the  THTTF  Itself  was  provided  in 
Section  2.  Discussions  of  details  of  the  instrumentation  systems, 
calibration  procedures  and  uncertainties  associated  with  the 
measurements  have  been  documented  in  Appendices  I-III. 

The  smooth  wall  fluid  dynamics  data,  presented  in  Section  3, 
showed  that  the  facility  was  operating  within  the  expected  bounds 
over  the  Reynolds  number  range  anticipated  in  future  rough  surface 
testing.  Skin  friction  coefficient  determinations  from  both  the 
Preston  tube  method  and  the  hot-wire  method  showed  excellent 
agreement  with  the  Schul tz-Grunow  correlation,  which  is  generally 
accepted  as  representing  smooth  wall  turbulent  skin  friction  data 
very  well.  The  agreement  shown  by  the  skin  friction  coefficients 
determined  using  the  hot-wire  method  was  especially  Important, 
since  that  is  the  technique  which  is  used  in  rough  wall  tests. 

Also,  the  behavior  of  the  hot-wire  data  was  consistent  with  the 
10-12%  estimate  of  the  uncertainty  in  skin  friction  coefficients 
determined  using  that  method. 

The  baseline  smooth  wall  heat  transfer  data,  presented  in 
Section  4,  showed  that  the  qualification  of  the  THTTF  was  completed 
and  successful  for  x-Reynolds  numbers  up  to  10  million.  Stanton 
numbers  determined  using  the  energy  balance  technique  were  In 
excellent  agreement  with  previous  smooth  wall  data  for  a  constant 
wall  temperature  boundary  condition,  which  had  been  reported  up  to 
x-Reynolds  numbers  of  3.5  million.  To  the  authors'  knowledge,  no 
smooth  wall  Stanton  number  data  has  been  reported  previously  in  the 
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3.5-10  million  range  covered  by  the  THTTF  data.  Stanton  numbers 
are  determined  In  the  THTTF  to  uncertainties  of  2-5%,  depending  on 
flow  conditions. 

Presented  In  Section  5  were  smooth  wall  Stanton  number  data 
taken  for  other  thermal  boundary  conditions:  (1)  constant  heat 
flux,  (2)  step  wall  temperature,  (3)  step  wall  heat  flux,  (4) 
linear  variations  in  wall  temperature,  and  (5)  bi-1 inear  variations 
in  wall  temperature.  These  data  were  the  first  to  be  reported,  to 
the  authors'  knowledge.  In  such  an  extended  x-Reynolds  number 
range. 

Presented  in  Section  6  were  the  Initial  rough  wall  data  taken 
using  the  THTTF.  These  data  were  for  a  surface  roughened  with 
hemispherical  elements  spaced  four  base  diameters  apart.  The 
Stanton  number  and  skin  friction  coefficient  data  were  compared 
with  predictions  from  the  discrete  element  method  using  unmodified 
the  roughness  models  developed  previous  to  this  effort.  The  skin 
friction  predictions  were  in  excellent  agreement  with  the  data. 

This  was  not  surprising  considering  the  large  size  of  the  data  base 
on  which  the  roughness  element  drag  coefficient  model  had  been 
developed.  The  heat  transfer  predictions  were  in  reasonably  good 
agreement  with  the  Stanton  number  data,  although  improvement  would 
be  desirable.  Since  the  roughness  element  Nusselt  number  model  was 
developed  based  on  heat  transfer  data  from  the  single  rough  surface 
used  in  the  Stanford  studies,  it  is  encouraging  that  the  heat 
transfer  predictions  are  as  good  as  they  are.  Further  development 
of  the  roughness  element  Nusselt  number  model  will  be  made  when 
data  are  taken  on  additional  rough  surfaces  in  the  THTTF. 
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APPENDIX  I 


EXPERIMENTAL  STANTON  NUMBER  DETERMINATION 
AND  ITS  UNCERTAINTY  ANALYSIS 

1.1  STANTON  NUMBER  DATA  REDUCTION  EQUATION 

The  Stanton  number  is  the  nondlmensional  convective  heat 
transfer  coefficient  and  may  be  defined  as 


where 

h  Is  the  convective  heat  transfer  coefficient 
p  is  the  density  of  freestream  air 

Cp  is  the  specific  heat  of  freestream  air 
U„  is  the  velocity  of  freestream  air 
The  rate  of  convective  heat  transfer  (q)  from  a  test  plate  to  the 
air  in  the  tunnel  is  defined  as 

q  -  hA(Tp  -  T0)  (1.2) 

where 

A  is  plate  area 

Tp  is  plate  surface  temperature 

Tq  is  freestream  air  total  temperature 

Solving  equation  (1.2)  for  the  convective  heat  transfer  coefficient 
and  substituting  Into  equation  (1.1)  gives  the  Stanton  number  at 
each  plate  as 


St 


q _ 

pCpU^(Tp  -  T0) 


(1.3) 
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In  order  to  determine  the  convective  heat  transfer  rate  q,  for 
each  plate,  the  corresponding  radiation  and  conduction  heat  losses 
are  required.  The  modes  of  heat  exchange  from  each  plate  are 
depicted  graphically  In  Figure  1.1.  Application  of  an  energy 
balance  to  a  plate  gives 

W  -  q  +  qc  +  qr  (1.4) 


where 

W  is  power  supplied  to  the  plate 
qr  is  radiation  heat  loss  rate 
qc  is  conductive  heat  loss  rate 
Solving  equation  (1.4)  for  the  convective  heat  transfer  rate  and 
substituting  into  equation  (1.3)  gives 


5t  W  -  qc  -  If 
PCpU*A(Tp  -  T0) 

The  radiation  heat  loss  rate  Is  modeled  using 
qr  -  oeA(Tp4  -  Tf4) 

where 

o  is  the  Stefan-Boltzmann  constant 
e  is  the  emissivlty  of  the  plate  surface 

Tr  is  the  freestream  recovery  temperature 

The  conduction  heat  loss  rate  Is  modeled  using 

qc  -  (UA)eff  (Tp  -  Trai1.) 


(1.5) 


(1.6) 


(1.7) 


where 

(UA)eff  Is  an  experimentally-determined  effective  conductance 
between  a  test  plate  and  the  side  rails 
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Trail  1s  the  side  rail  temperature  at  the  axial  location  of 
the  plate 

Substitution  of  (1.6)  and  (1.7)  into  (1.5)  gives  the  final 
form  of  the  data  reduction  expression  for  the  Stanton  Number 


W  -  (UA)eff  (Tp  -  Trai])  -  oeA(Tp  -  Tr*) 
pCpU_A(Tp  -  T0) 


(1.8) 


This  expression  shows  explicitly  most  of  the  variables  in¬ 
volved  In  the  experimental  Stanton  number  determination.  Addi¬ 
tional  variables  enter  In  the  determination  of  the  static  and  total 
temperature  of  the  freestream  air  and  In  the  moist  air  property 
calculations  for  Cp  and  p.  The  freestream  air  total  and  static 
temperatures  are  calculated  using  the  measured  recovery  temperature 
and  a  recovery  factor,  R,  for  the  probe: 


(1.9) 

(1. 10) 


The  functional  relationship  for  the  moist  air  specific  heat  calcu¬ 
lation  is 

T 

Cp  -  Cp(T„,  Twb,  Pbar,  Cpai(_,  CpH^0) 


where 

T_  is  the  freestream  air  static  temperature  and  is  also 
taken  as  the  dry  bulb  temperature 
Twb  is  the  freestream  air  wet-bulb  temperature 
Pbar  Is  the  barometric  pressure 

Cpair  the  dry  air  spec*f*c  heat 
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cPh20 


Is  the  water  vapor  specific  heat 


The  functional  relationship  for  the  moist  air  density  is 


P  “  P(T„.  Twjj,  Pbar) 

The  Stanton  number  determination  for  each  plate,  therefore, 
Involves  the  following  thirteen  variables  which  are  either  measured 
or  found  from  a  reference  source: 

•  Plate  heater  power  (W) 

•  Recovery  temperature  (Tr) 

•  Plate  temperature  (Tp) 

•  Rail  temperature  (Tra^) 

•  Wet-bulb  temperature  (Twb) 

•  Effective  conductance  ((UA)eff) 

•  Plate  area  (A) 

•  Barometric  pressure  (Pbar) 

•  Specific  heat  of  dry  air  (cpa1r) 

•  Specific  heat  of  water  vapor  (CpH  Q) 

•  Freestream  air  velocity  (U„)  2 

•  Recovery  factor  (R) 

•  Emissivity  (c) 


1.2  UNCERTAINTY  ANALYSIS  OVERVIEW 

The  detailed  uncertainty  analysis  procedure  follows  Coleman 
and  Steele  1989  [37],  which  Is  consistent  with  the  1986  ANSI/ASME 
Standard  on  Measurement  Uncertainty  [52].  Briefly,  the  true  value 
of  a  quantity,  which  Is  approximated  by  an  experimental  result  r, 
lies  within  the  Interval  r  ±  Uf  with  95%  confidence.  Here  Ur  Is 
the  uncertainty  In  the  result  determined  from  the  root-sum-square 
combination  of  the  bias  limit  of  the  result,  Br,  and  the  precision 
limit  of  the  result,  Pr, 

Ur  -  (Br2  +  Pr2)1/2 
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(I. 11) 


For  a  result  (such  as  St)  which  Is  a  function  of  J  variables 
and  parameters 

r  -  r(Xj,  X2 . Xj)  (1.12) 

the  propagation  of  the  precision  limits  Px.  of  the  measured 
variables  into  the  result  Is  given  by 

>v-n  £-px,)2]I/2 

1-1  axl  1 

and  the  propagation  of  the  bias  limits  of  the  variables  Into 
the  result  Is  given  by 


As  explained  In  [37],  there  Is  a  term  such  as  the  second  one  on  the 
right  hand  side  of  equation  (1.14)  for  each  pair  of  variables  that 
have  portions  (B1^  and  B1^)  of  their  bias  limits  which  are 
perfectly  correlated.  These  usually  arise  when  the  transducers 
used  to  measure  different  variables  have  been  calibrated  against 
the  same  standard  or  when  two  variables  (often  temperatures  or 
pressures)  are  measured  with  the  same  transducer. 

In  the  current  experiments,  all  of  the  Px^s  are  negligible 
compared  to  the  bias  limits,  so 

Pst  -  0  .  (1.15) 

Exceptions  to  this  occur  at  low  freestream  velocities— U„  i  12 
m/sec  for  the  smooth  wall  tests  and  U„  i  6  m/sec  for  the  rough  wall 
tests— for  which  the  heat  transfer  coefficients  are  relatively  low. 
At  these  conditions  the  time  constant  of  the  THTTF  Is  large  enough 
so  that  the  relatively  long  period  variations  In  facility  line 


112 


voltage  to  the  test  plate  heater  circuits  and  In  the  temperature  of 
the  Incoming  make-up  water  for  the  heat  exchanger  loop  affect  the 
ability  to  hold  a  tight  steady  state  condition.  These  annoyances 
could  be  overcome  with  additional  expenditures  for  power  condition¬ 
ing  equipment  and  a  water  chiller  system;  however,  the  observed 
run-to-run  scatter  in  St  results  at  these  low  velocities  is  within 
acceptable  limits.  Observations  of  the  St  results  for  eight  U„  - 
12  m/sec  replications  with  the  smooth  wall  and  three  U.  -  6  m/sec 
replications  with  the  rough  wall  produced  a  95%  confidence  estimate 
of  P5t  -  3%  for  those  conditions.  This  is  present  because  of 
system  unsteadiness  and  not  because  of  measurement  uncertainty. 

Application  of  equation  (1.14)  to  the  case  of  equation  (1.8) 

gives 


f 3St,o  o  f3St.?  o 
- ZB  t  +  - ZB  t  +  , 

l3TpJ  TP  l3TrJ  Tr  3Trai 1 


(  8St  l2P2 


]*bzt  +  (— ]VT 

'rail  l3TwbJ  ' 


3St 


9Pbar 

JCB‘ 

asti2, 

r2 

3A  1 

8  A 

.  3St 

f3St 
’bar  +  3W  ^ 


2o  2 


-3St, 


2d2 


r9St, 


3C, 


3St  . j  o 

-)ZBZr_  ♦  t 


3(UA)eff 


pa1r  air 

]2B2(UA)eff  *  g)2*' 


2d2 


8St  l2B2 


3C 


wb 

R 

(1.16) 


pH20 


h2o 


f 3St. f3St,  . 3St, .  3St  , 

+  2  -  - B't  B't  +  2 -  - B't  B' 

l3T  Jl3T  '  n  lr  1»t  ilaT  _ J  in 


'p  ■  r 


V  Tra11 


f 3St, .  3St  . 

+  2  -  - )B't  B't 

9Tr  9Tra11  r  r»" 


where,  as  discussed  later,  the  only  portions  of  biases  that  are 
considered  correlated  are  those  arising  from  calibrating  the  therm¬ 
istors  which  measure  Tp,  Tr,  and  Tra11  against  the  same  reference 
standard. 
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1.3  MEASUREMENTS  AND  UNCERTAINTIES 

The  experimental  determination  of  Stanton  number  requires 
values  for  thirteen  variables,  some  of  which  can  be  measured  di¬ 
rectly  and  some  which  cannot.  Therefore,  the  methodology  used  In 
determination  of  each  variable  will  be  discussed.  Discussion  of 
the  determination  of  each  variable  will  include  description  of  the 
required  measurement  system,  Information  about  the  calibration 
procedures,  and  uncertainty  estimates  for  each  variable. 

1.3.1  Plate  Heater  Power 

The  power  supplied  to  the  heater  pad  of  each  plate  Is  measured 
by  a  high  precision  ac  watt  transducer  coupled  to  a  Hewlett  Packard 
(HP)  Model  3054  A  Automated  Data  Acquisition  and  Control  System 
(ADACS),  which  in  turn  Is  connected  to  a  Model  220  microcomputer. 
This  watt  transducer  Is  a  single  phase  transducer  with  a  rated 
output  of  1  ma  corresponding  to  500  watts.  The  manufacturer  speci¬ 
fies  ±  0.2%  of  reading  accuracy  and  0  to  1  ma  dc  current  output 
proportional  to  electrical  power. 

A  separate  power  circuit  is  used  for  each  Individual  plate 
heater.  There  are  24  Identical  power  circuits  for  the  24  test 
plates.  A  single  watt  transducer  (Ohio  Semitronics  Inc.  Model 
EW5-B)  is  used  for  all  power  measurements.  The  power  delivered  to 
a  plate  heater  is  measured  by  routing  the  power  through  the  watt 
transducer  by  switch  closures  using  the  ADACS.  Since  the  ADACS 
cannot  process  current  signals  directly,  the  transducer's  output  is 
measured  indirectly.  A  7.5  Kfl  resistor  is  shunted  across  the 
transducer's  output  lines  so  that  the  current  output  is  transformed 
into  a  measurable  voltage.  The  shunt  resistor  Is  sized  to  compen¬ 
sate  for  the  small  current  output  from  the  watt  transducer.  The 
current  output  (ma)  from  the  watt  transducer  is  then  obtained  using 
Ohm's  law 


1 


(-]  (J_ 

lR  1 1000 


(1.17) 
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This  current  output  Is  then  translated  Into  power  (watts)  using  the 
relationship 

W  -  500  x  1  (1.18) 

A  calibration  plate  heater  circuit  was  used  to  check  the  cali¬ 
bration  of  the  watt  transducer  by  comparison  of  the  transducer 
measurement  to  the  heater  power  (Wact)  determined  using  the  ADACS. 
This  power  was  obtained  by  measuring  the  ac  voltage  drop  across  the 
plate  heater,  the  resistance  of  the  plate  heater,  and  using 


Due  to  the  importance  of  resistance  and  voltage  measurement  in 
determination  of  both  transducer  and  actual  powers,  extreme  care 
was  exercised  to  utilize  the  ADACS  properly.  In  particular,  the 
resistances  from  the  plate  heater  and  shunt  resistor  were  measured 
using  the  four-wire  technique.  In  this  method,  the  resistance  of 
the  transmission  line  Is  measured  and  Is  subtracted  from  the  meas¬ 
ured  total  resistance  automatically.  Therefore,  the  resistance 
obtained  by  the  four-wire  technique  represents  the  load  resistance 
alone. 

Uncertainties:  The  high  resolution  digital  voltmeter  (3456A) 
used  for  both  voltage  and  re  *  ~e  measurements  has  a  voltage 
accuracy  of  ±  0.007%  of  reading  with  an  ADACS  environment  tempera¬ 
ture  of  23  C  and  an  additional  0.0002%  error  for  every  1  C  In  the 
environment  temperature  above  or  below  23  C.  Resistance  accuracy 
for  the  ADACS  Is  0.008%  of  reading,  and  the  environment  temperature 
coefficient  Is  0.004%/C  for  four-wire  ohm.  It  was  assumed  that  the 
manufacturer's  specification  on  the  accuracies  of  the  resistance 
and  voltage  measurements  are  the  bias  limits  with  95%  confidence 
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level  so  that  averaging  multiple  readings  would  not  reduce  these 
estimates.  Therefore,  the  bias  limits  on  the  voltage  and  resis¬ 
tance  measurements  become 

—  -  ±  7  x  10"5 

V 

Bn  c 

—  -  ±  8  x  icrb 

R 

The  bias  limit  In  the  determination  of  the  watt  transducer 
output  current  using  equations  (1.14)  and  (1.17)  is 


ft2 


±  0.01% 


and  the  bias  limit  In  the  determination  of  the  actual  power  mea¬ 
surement  using  equations  (1.14)  and  (1.19)  is 


(Vt]2  .  (!r)2  +  (2  !v)2 

“act  R  V 

By 

(_££!)  .  ±  o.02% 

"act 

These  are  so  small  as  to  be  negligible  for  our  purposes.  Precision 
errors  were  also  observed  to  be  negligible. 

Although  the  output  from  the  watt  transducer  can  be  measured 
with  high  accuracy  as  indicated  above,  how  well  this  Is  translated 
Into  power  using  equation  (1.18)  must  be  determined  by  comparing 
the  power  Indicated  by  the  watt  transducer  to  the  actual  power  as 


116 


found  from  the  calibration  tests.  (The  Inductance  In  the  heating 
element  was  accounted  for,  with  the  power  factor  being  greater  than 
0.999  [23].)  The  power  Indicated  by  equation  (1.18)  from  the  output 
of  the  watt  transducer  was  compared  with  the  actual  power  as  meas¬ 
ured  by  the  ADACS  (1.19)  using  172  points  over  the  0-250  watt  range 
of  Interest.  Figure  1.2  shows  the  cumulative  probability  curve  of 
the  absolute  values  of  the  percent  differences.  As  shown,  a  95% 
confidence  estimate  of  the  uncertainty  In  the  watt  transducer  power 
measurement  based  on  these  points  Is  ±  0.9%  of  reading.  This 
appears  as  a  bias  error  In  W  when  equation  (1.18)  Is  used. 

1.3.2  Temperatures 

Temperatures  are  measured  using  thermistors,  which  are  tem¬ 
perature  sensitive  resistors  with  a  negative  temperature  coeffi¬ 
cient.  These  thermistors  have  a  nominal  resistance  of  50,000  ohms 
at  25  C  and  are  highly  sensitive  to  small  temperature  changes 
(about  1-2  Kfl/C).  They  are  guaranteed,  by  the  manufacturer,  to 
have  ±  0.2  C  Interchangeability  over  a  range  of  temperatures  from  0 
C  to  70  C.  The  resistances  of  the  thermistors  are  measured  by  the 
ADACS.  These  thermistors  are  used  to  determine  the  freestream  air 
temperature,  the  test  plate  temperatures,  and  the  metal  support 
rail  temperatures. 

The  calibration  of  the  thermistors  was  done  In  a  Blue  M  Model 
MR-3210A-1  constant  temperature  bath.  The  bath  temperature  was 
monitored  by  a  Hewlett  Packard  Model  2804A  quartz  thermometer 
Instrumented  with  a  Model  1811 1A  quartz  probe.  The  absolute  accu¬ 
racy  of  this  thermometer/probe  combination  Is  specified  by  the 
manufacturer  as  l  0.040  C  over  a  range  of  -50  C  to  150  C.  The 
thermistors  were  placed  Individually  Inside  glass  test  tubes  to 
protect  them  and  avoid  their  contamination.  To  ensure  effective 
conduction  of  heat  from  the  water  bath  to  each  thermistor,  each 
test  tube  was  filled  with  Mega therm  201  (by  Omega  Engineering, 
Inc.),  which  Is  a  high  thermal  conductivity,  filled  silicone  paste. 
The  spatial  variation  In  the  temperature  of  the  bath  was  found  to 
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be  about  ±  0.4  C.  This  variation  was  minimized  to  l  0.02  C  by 
centering  the  test  tubes  containing  thermistors  around  the  quartz 
probe  In  groups  of  fourteen.  Since  the  reliability  of  measured 
thermistor  resistances  depends  on  the  accuracy  of  measurements  made 
by  the  ADACS,  the  proper  use  of  the  ADACS  during  the  measurement 
process  was  carefully  examined. 

The  thermistor  calibrations  were  performed  for  the  range  of 
temperatures  22  C  -  50  C  using  six  evenly  spaced  points  over  this 
range.  The  thermistors  are  extremely  nonlinear  but  their  behavior 
can  be  very  closely  approximated  by  the  Stelnhart-Hart  equation  as 

T[K]  -  1 / [A  +  B  in  R  +  C(tn  R)3]  (1.20) 

where  R  Is  resistance  In  ohms.  The  curveflt  constants  A,  B  and  C 
were  calculated  using  the  thermistor  manufacturer's  data  as  A  > 
9.6401  x  10-4,  B  -  2.1095  x  10-4  and  C  -  8.48  x  10'8. 

Uncertainties:  The  temperatures  obtained  using  the  measured 
thermistor  resistances  In  the  Stelnhart-Hart  relation  (equation 
(1.20))  were  compared  with  the  temperatures  from  the  quartz  ther¬ 
mometer.  The  difference  between  the  temperatures  measured  by  the 
quartz  thermometer  and  the  temperatures  calculated  from  the  meas¬ 
ured  thermistor  resistances  using  the  curvefit  equation  was  deter¬ 
mined  for  a  total  of  360  calibration  points.  Figure  1.3  shows  the 
cumulative  probability  curve  of  the  absolute  values  of  the  differ¬ 
ences.  As  shown,  95%  are  less  than  about  0.09  C.  This  Is  taken  as 
a  bias  limit  since  precision  errors  were  observed  to  be  negligible. 

The  bias  limits  for  the  elemental  error  sources  which  affect 
the  various  thermistors  are  estimated  as 


Elemental  error  source  Bias  limit 
Calibration:  Quartz  probe  ±  0.04  C 
Calibration:  Bath  nonunlformlty  &  curveflt  ±  0.09  C 
Installation  In  test  plates  ±  0.1  C 
Nonuniformity  In  side  rail  temperatures  l  0.4  C 
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The  Installation  error  for  the  test  plate  thermistors  was 
estimated  based  on  the  temperature  variations  predicted  within  a 
test  plate  during  the  design  calculations  by  Norton  [22].  Since 
the  temperature  Indicated  by  the  plate  thermistors  Is  used  as  the 
plate  surface  temperature,  Tp,  In  equation  (1.8),  such  variations 
are  the  source  of  a  bias  error.  The  overall  bias  limit  for  each 
plate  temperature  measurement  Is  found  from  the  root-sum-square  of 
the  three  appropriate  elemental  error  sources  as 

m 

BTp-  [(0.04)2  +  (0.09)2  +  (0.1)2]172  «  0.14  C 

For  measurements  of  the  side  rail  temperature,  Trail,  at  a 
given  axial  location,  the  overall  bias  limit  Is  found  from  the 
root-sum-square  of  the  three  appropriate  elemental  error  sources  as 

BTrail  -  [(0.04)2  +  (0.09)2  +  (0.4)2]172  -  0.4  C 

For  measurements  of  the  freestream  air  recovery  temperature, 
Tr,  no  Installation  bias  appears  and  so 

Byr-  [(0.04)2  +  (0.09)2]172  -  0.1  C 

Since  the  same  calibration  standard  (quartz  probe)  was  used 
for  all  thermistors,  the  standard's  bias  limit  Is  a  correlated  bias 
error  source  that  has  to  be  accounted  for  in  the  data  reduction 
uncertainty  analysis.  Thus,  B'y^  •  B'yr  -  B'Tra^  -  0.04  C  in 
equation  (1.16). 

1.3.3  Effective  Conductance 

The  effective  conductance  (UA)eyy  was  determined  from 
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(UA)eff 


(1.21) 


(Tp  -  Trai1> 


where 

(UA)eff  U  the  effective  conductance  between  a  test  plate  and 
the  support  rails 

qc  Is  the  Input  power  to  a  plate  heater  under  conduct¬ 

ance  calibration  conditions 

In  order  to  estimate  the  effective  conductance,  an  experimental 
approach  was  taken.  Insulation  was  placed  over  the  top  of  the  test 
plates,  which  then  were  heated  by  power  Input  to  the  plate  heaters. 
Since  there  were  no  radiative  or  convective  heat  losses  from  the 
covered  plates,  the  total  Input  power  to  the  heater  plates  was 
equal  to  the  conductive  heat  transfer  loss  from  the  plates,  which 
was  modeled  with  equation  (1.21). 

The  support  rails  were  heated  by  two  150  watt  tape  heaters  In¬ 
stalled  on  each  side  rail.  When  the  plate  temperatures  and  the 
rail  temperatures  reached  a  prespecified  temperature,  the  powers  to 
the  plate  heaters  were  reduced  to  0.6  watts  and  the  rail  heaters 
were  turned  off.  The  test  plate  temperatures  and  the  rail  tempera¬ 
tures  at  x  -  1.2  m  were  monitored  until  the  temperature  differences 
between  each  test  plate  and  the  rail  temperature  approached  a 
steady  state  condition. 

For  each  plate,  from  each  temperature  difference  and  the  cor¬ 
responding  Input  power,  an  effective  conductance  was  determined.  A 
single  value  of  effective  conductance  (0.42  watt/C)  was  determined 
for  use  with  all  test  plates.  This  value  was  obtained  by  averaging 
the  effective  conductances  of  test  plates  6  through  18.  The  preci¬ 
sion  Index  S  of  this  sample  of  effective  conductances  was  0.08 
watt/C  . 

Uncertainties:  The  uncertainties  were  determined  by  using  the 
uncertainty  analysis  equations  (1. 11)  -  (1.14)  and  equation  (1.21) 
which  models  the  conductive  heat  loss.  Since 
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(u*)eff 


qc 

Tp  ”  Tra11 

then 


(UA) ef f j2  „  fBqcj2  +  r  Btp  ,]2  +  r  Byra1 1  y 
(UA)eff  qc  (Tp  -  Tra^)  (Tp  -  Trail) 

2  B't  B*t  ( - - - ]( - - - 

p  rail  Tp  -  Tra^  Tp  -  Tpa1 


From  previous  discussion  we  have 


Bt  «  0.14  C 
TP 


B' 


B'Tra1, 


0.04  C 


and 


qc 

—  -  0.009 
qc 

The  bias  limit  for  Tra11  does  not  Include  the  0.4  C  nonuniformity 
contribution  since  only  the  middle  13  plates  were  used  In  the 
calibration  test  to  avoid  the  larger  variations  In  Tra11  at  the 
beginning  and  end  of  the  test  section.  For  this  case,  then, 

8Tran  *  "-I  ‘ 

Using  a  nominal  value  for  (Tp  -  Tra11)  of  1.3  C  from  the  calibra¬ 
tion  test,  the  bias  limit  In  the  effective  conductance  becomes 
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B(UA)effl2 
(UA>eff  J 


(0.009)2  +  [j£]2  *  (y^)2  -  2(0.04) (0.04) 


B(UA)eff 

(UA)eff 


]  -  0.13 


Using  the  precision  limit  (P  -  tS  -  0.174  watt/C)  at  a  95% 
confidence  level  calculated  from  the  t-distribution  (n-13,  t-2.179) 
and  the  precision  index  S,  the  overall  uncertainty  for  the  effec¬ 
tive  conductance  becomes 


fU(UA)effl2 

(UA)eff  J 

fU(UA)eff1? 
l(UA)eff  J 


B(UA)eff  2  P(UA)eff  2 

‘tsawt!2  +  lm^r] 

(0*13)2  +  *  (0*13)2  +  (0*42)2 


fU(UA)eff 
(UA)ef f 


]  -  45% 


Since  the  constant  value  of  0.42  watt/C  is  always  used  for  (UA)eff, 
this  uncertainty  Is  fossilized  [37]  into  a  bias  limit  when  (UA)ef^ 
is  used  in  the  calculation  of  Stanton  numbers.  Therefore, 

B(UA)eff  "  45%* 


1.3.4  Area 

The  surface  area  of  the  plates  Is  determined  from  the  lengths 
of  the  sides  (10.16  cm  by  45.72  cm) 

A  -  (Sj ) (S2) 

The  plates  were  manufactured  with  length  and  width  tolerances  of 
±  0.0025  cm. 
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Uncertainties:  The  uncertainty  In  the  plate  area  Is  assumed 
to  be  all  bias  and  may  be  expressed  as 


Substitution  of  the  plate  dimensions  and  bias  limits  gives 


0.0025,7  ,0.0025 

- K  +  I - 

10.16  J  l45.72 


2 


,B» 

(— ]  -  0.03* 


The  possibility  of  thermal  expansion  of  the  test  plates  (which 
could  introduce  additional  uncertainty  into  the  area)  was  consid¬ 
ered.  It  was  found  that  this  effect  was  negligible  compared  to  the 
bias  due  to  manufacturing  tolerance. 

1.3.5  Air  Density  and  Specific  Heat 

The  fluid  in  the  test  section  is  actually  a  mixture  of  dry  air 
and  water  vapor.  Therefore,  fluid  properties  such  as  density  and 
specific  heat  for  the  test  air  will  depend  on  the  ratio  of  dry  air 
and  water  vapor  In  the  mixture.  The  ratio  of  dry  air  and  water 
vapor  in  air  is  reflected  by  the  partial  pressures  of  each.  The 
density  and  specific  heat  of  the  mixture  may  be  expressed  in  terms 
of  the  partial  pressures  of  the  dry  air  and  water  vapor.  Once  the 
partial  pressures  of  the  dry  air  and  water  vapor  are  known,  they 
may  be  substituted  Into  expressions  for  the  density  and  specific 
heat  of  the  mixture. 

Psychrometrlcs:  The  partial  pressure  of  the  water  vapor  at 
saturation,  evaluated  at  the  dry  bulb  temperature  T„,  Is  given  by 


123 


•  Pws  *  exp[C8/T.  +  C9  +  (C10)T.  +  (CU)T.2  +  (C12)T„3 
+  (C13)  ln(TJ] 

and  the  partial  pressure  of  the  water  vapor  at  saturation,  evalu¬ 
ated  at  the  wet  bulb  temperature,  is  given  by 

PWswb  *  exp[C8/Twb  +  C9  +  (C10)Twb  +  (Cll)Twb2  +  (C12)Twb3 
+  (C13)  *n(Twb)] 

where  the  constants  are 

C8  -  -10440.4 
C9  -  -11.2946669 
CIO  -  -0.02700133 
Cll  -  1.289706  x  10“5 
C12  -  -2.478068  x  10~9 
C13  -  6.5459673 


the  temperatures  are  In  degrees  R,  and  the  partial  pressures  in 
psia  [53].  Pwswb  can  be  used  to  determine  the  humidity  ratio  at 
wet  bulb  temperature  from 


Wwb 


(0.62198)  Pwswb 
(pbar  "  pwswb) 


The  humidity  ratio  is  then  obtained  as 

[1093  -  (0.556)1^30^)  -  0.24(T.  -  Twb) 
1093  +  (0.44)T.  -  Twb 

and  the  humidity  ratio  at  saturation  as 

(0.62198)  Pws 

u  ■  - 

pbar  ~  pws 
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The  degree  of  saturation  Is  defined  as 


and  the  relative  humidity  Is  then  calculated  using 

_ Us _ 

*  '  i  -  (i  -  UslOWW) 

The  partial  pressure  of  the  water  vapor  Is  determined  from  the 
relative  humidity  and  the  partial  pressure  of  the  water  vapor  at 
saturation,  thus 


pw  *  pws  ^ 

Once  the  partial  pressure  of  the  water  vapor  and  the  ambient 
pressure  are  known,  Dalton's  law  of  partial  pressure  may  be  applied 
to  determine  the  partial  pressure  of  the  dry  air  as 

pa  *  pbar  "  pw 

Density:  The  density  of  the  moist  air  flowing  through  the 
test  section  may  be  expressed  as 


ma  +  mw 


'm 


m: 


a  |  ^ 


fm 


'm 


The  Ideal  gas  law  may  be  applied  to  the  dry  air  and  water  vapor  and 
written  as 


P  V 
ravm 

P  V 
rwTm 


naRT- 

nwRT. 
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The  number  of  moles  of  a  gas  is  equal  to  the  quotient  of  the  mass 
of  gas  divided  by  the  molecular  weight  of  the  gas.  Therefore,  the 
number  of  moles  of  dry  air  and  water  vapor  may  be  expressed  as 


nw  * 

The  specific  gas  constants  of  dry  air  and  water  vapor  are  defined 
as 


Ra  -  R/Ma 
R*  ■  R/N* 

Substitution  into  the  Ideal  gas  law  relations  gives 


pavm 

pwvm 


maRaT- 

W- 


Solving  for  the  dry  air  mass  and  the  water  vapor  mass  and 
substituting  in  the  density  equation  yields 


P 


(1.22) 


This  equation  states  that  the  density  of  a  mixture  of  gases  Is  the 
sum  of  the  densities  of  the  component  gases  at  their  respective 
partial  pressures. 

Specific  Heat:  The  specific  heat  of  the  test  air  may  be  de¬ 
termined  by  applying  an  energy  balance  to  the  dry  air  and  water 
vapor  system. 


cp%T-  "  cpamaT-  +  Cpv/V- 
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Solving  for  the  specific  heat  of  the  mixture  and  substituting  for 
the  mass  of  mixture  yields 


cpa  ma  +  cpw  "V 

C  .  - 

P  ma  +  % 

Substituting  for  the  mass  of  the  dry  air  and  the  mass  of  water 
vapor  and  dividing  by  the  total  number  of  moles  In  the  mixture 
gives 

cPa  Ma"a/n  cPw  W" 

CP  "  Ma(r»a/n)  +  Mw(nw/n)  +  Ma(na/n)  +  Mw(nw/n) 

From  the  thermodynamic  principles  of  the  mixture  of  gases,  we  know 
that  the  ratio  of  each  partial  pressure  to  the  total  pressure  is 
equal  to  the  mole  fraction  of  each  component.  Therefore,  we  have 


P  n 


Substitution  for  the  number  of  moles  In  terms  of  partial 
pressures  gives  the  final  expression  for  the  specific  heat  of  the 
test  air  as 


C 


P 


Cpa  Va  +  ‘p*  V* 

MaPa  ♦  Vw 


(1.23) 


It  should  be  noted  that  the  following  assumptions  have  been 
made  In  this  section.  The  curveflt  equations  used  to  obtain  par¬ 
tial  pressures  are  assumed  very  accurate,  and  their  uncertainties 
are  negligible  In  comparison  with  the  measurement  uncertainties. 
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The  ideal  gas  thermodynamic  relationships  used  to  calculate  the 
test  air  properties  are  applicable  to  the  dry  air  and  saturated 
water  vapor,  and  the  uncertainties  Introduced  by  their  application 
to  a  non-ideal  gas  are  negligible  as  compared  with  the  measurement 
uncertainties  and  the  uncertainties  In  the  properties  obtained  from 
reference  tables. 

Uncertainties:  A  sling  psychrometer  is  inserted  through  an 
access  hole  in  the  top  of  the  test  section  and  is  used  to  find  the 
wet  bulb  temperature  of  the  test  air.  The  bias  limit  in  Twb  is 
taken  as  1.0  C  and  precision  error  is  assumed  negligible.  The 
static  air  temperature  determined  from  the  recovery  temperature 
measured  by  the  thermistor  probe  Is  used  for  the  dry  bulb  tempera¬ 
ture.  A  temperature  compensated  barometer  with  a  resolution  of  0.5 
mm  Hg  is  used  to  determine  the  ambient  pressure.  Since  the  pres¬ 
sure  of  the  test  air  mixture  (dry  air  +  water  vapor)  in  the  test 
section  is  kept  equal  to  the  outside  ambient  pressure,  using 
make-up  air  controls,  the  pressure  read  from  the  barometer  is  used 


as  the  air  mixture  pressure.  A  bias  limit  of  1.0  mm  Hg  and  negli¬ 
gible  precision  error  are  assumed  for  Pbar. 

The  specific  heats  for  dry  air  and  for  water  vapor  were  ob¬ 


tained  as  [54]  Cpa1r  -  1.006  kJ/kg  C  and  C 
an  uncertainty  in  each  estimated  as  0.5%. 


d  *  1.86  kJ/kg  C  with 

“HpU 

These  uncertainty  esti¬ 


mates  are  fossilized  into  bias  limits  when  the  specific  heat 


values  are  used  in  calculations. 


1.3.6  Air  Velocity 

A  Pitot  probe  Is  used  to  determine  the  freestream  velocity. 
The  Pitot  probe  channels  the  stagnation  pressure  and  static  pres¬ 
sure  exerted  by  the  freestream  Into  a  differential  pressure  trans¬ 


ducer  so  that  the  freestream  dynamic  pressure  can  be  measured. 
Once  the  dynamic  pressure  (AP)  Is  known,  the  freestream  velocity 
may  be  determined  from 
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(1.24) 


U.  -  l^)1'2 

where  p  Is  the  density  of  the  freestream  air.  The  procedure  for 
determination  of  the  test  air  density  has  already  been  described. 
The  dynamic  pressure  determination  methodology  Is  presented  next. 

Dynamic  pressure  measurement  is  performed  using  two  differen¬ 
tial  pressure  transducers  with  ranges  of  0.55  and  3.45  KPa.  These 
transducers  cover  the  full  range  of  dynamic  pressures  expected. 
Their  accuracy,  specified  by  the  manufacturer,  is  ±  0.5  percent  of 
full  scale.  Each  transducer  provides  a  voltage  output  of  0-5  Vdc 
proportional  to  the  applied  pressure  difference.  The  voltage 
outputs  of  the  pressure  transducers  are  measured  by  the  ADACS. 

The  calibration  of  each  pressure  transducer  was  accomplished 
by  employing  a  very  sensitive  water  micromanometer  as  the  pressure 
source.  The  25.4  cm  range  micromanometer  Is  equipped  with  a 
magnifier  which  amplifies  the  fluid  meniscus  at  the  reference 
hairline  and  provides  direct  reading  Indication  to  0.0025  cm  of 
water. 

Each  pressure  transducer  was  calibrated  individually.  Various 
pressures  within  the  pressure  transducer  range  were  generated  using 
the  micromanometer  and  were  applied  to  the  transducer.  The  values 
of  these  generated  pressures  indicated  by  the  micromanometer  were 
recorded  and  the  corresponding  voltage  outputs  from  the  pressure 
transducer  were  measured  using  the  ADACS  in  a  manner  corresponding 
to  actual  testing.  The  number  of  pressure  calibration  points 
obtained  for  the  0.55  and  3.45  KPa  range  transducers  were  23  and 
10,  respectively.  Both  pressure  transducers  demonstrated  small  but 
stable  voltage  outputs  at  zero  pressures  (zero  shift).  The  pres¬ 
sure  calibration  data  collected  from  each  transducer  was  corrected 
for  the  zero  shift. 
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Subsequently,  the  data  of  each  transducer  was  used  to  arrive 
at  an  appropriate  curvef It  equation  for  that  transducer.  A  linear 
curvefit  equation  for  the  3.45  KPa  range  pressure  transducer  was 
satisfactory,  but  a  quadratic  equation  was  necessary  to  fit  the 
0.55  KPa  range  pressure  transducer  calibration  data  satisfactorily. 

Uncertainties:  A  comparison  between  the  direct  pressure  meas¬ 
urement  data  (micromanometer  readings)  and  the  pressures  calculated 
from  the  curvefit  equation  was  made  for  each  pressure  transducer. 
The  results  Indicated  that  the  bias  limits  associated  with  using 
the  calibration  curvefits  for  the  0.55  and  3.45  KPa  range 
transducers  were  0.5  percent  and  0.1  percent  of  reading,  respec¬ 
tively.  Precision  errors  were  observed  to  be  negligible  during  the 
calibration  process,  and  the  bias  error  Inherent  In  the  micro- 
manometer  (which  was  used  as  the  calibration  standard)  was  assumed 
negl Igible. 

During  the  actual  dynamic  pressure  measurement  when  a  Pitot 
probe  in  the  flow  stream  Is  the  pressure  source  for  the  transducer, 
the  elemental  bias  error  sources  are  the  Pitot  probe  and  the  meas¬ 
urement  system  calibration.  Therefore,  the  calibration  bias  and 
the  additional  biases  Introduced  by  the  Pitot  probe  must  be  com¬ 
bined  using  RSS  to  arrive  at  the  overall  bias  limit  for  actual 
dynamic  pressure  measurements.  Biases  in  the  dynamic  pressure  due 
to  errors  caused  by  Pitot  probe  design,  use,  and  misalignment  have 
been  estimated  at  0.5  percent  since  the  freestream  flow  is  uniform 
and  relatively  free  of  perturbations  and  since  the  Pitot  probe  Is 
very  carefully  aligned  with  the  flow.  The  bias  limit  estimates  are 

for  Pitot  probe  0.5%  of  reading 

for  0.55  KPa  transducer  (calibration)  0.5%  of  reading 
for  3.45  KPa  transducer  (calibration)  0.1%  of  reading. 

Therefore,  bias  limits  for  dynamic  pressure  measurements  become 


130 


0.7% 


(— )  -  [(0.005)2  +  (0.005)2]172  - 
AP 

for  the  0.55  KPa  transducer,  and 

g 

(— )  -  [ (0 . 001 ) 2  +  (0.005)2]1/2  -  0.5% 

AP 

for  the  3.45  KPa  transducer. 

The  bias  limit  in  the  free  stream  velocity  is  dependent  on  the 
bias  limit  in  the  dynamic  pressure  and  on  the  bias  limit  In  the 
density.  As  discussed  above,  bias  limits  for  the  dynamic  pressures 
measured  using  the  0.55  KPa  and  3.45  KPa  transducers  are  0.5%  and 
0.7%  of  reading,  respectively.  For  the  bias  limit  of  the  air 
density,  based  on  an  analytical  uncertainty  analysis  using  nominal 
values  of  the  related  variables,  an  estimated  value  of  0.3%  is 
used.  This  Is  a  conservative  value  and  does  not  change  signifi¬ 
cantly  with  small  changes  In  the  environment  and/or  operating 
conditions  of  the  tunnel. 

Applying  equation  (1.14)  to  equation  (1.24),  the  expression 
for  the  bias  limit  of  the  freestream  velocity  is 


which,  after  substitution  of  the  appropriate  partial  derivatives 
and  division  by  the  velocity  gives 


The  bias  limit  when  the  0.55  KPa  pressure  transducer  is  used  is 
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,eU. 

(jjH  -  0.4* 

The  bias  limit  when  the  3.45  KPa  pressure  transducer  is  used  is 
BU 

•ur)2  ■  °-3* 

Precision  errors  have  been  observed  to  be  negligible  In  the 
determination  of  U„,  so  the  precision  limit  for  U„  is  assumed  zero. 

1.3.7  Air  Temperature  Probe  Recovery  Factor 

The  temperature  obtained  from  the  air  thermistor  is  the  recov¬ 
ery  temperature  Tr  The  freestream  total  temperature  and  static 
temperature  are  computed  from  Tr  using  equations  (1.9)  and  (I. 10), 
but  a  value  for  recovery  factor  R  is  necessary  for  the  calculation. 
Based  on  the  review  of  available  data  [55],  a  value  of  R  -  0.86  was 
chosen  with  an  uncertainty  of  0.09.  This  uncertainty  is  fossilized 
Into  a  bias  limit  when  the  recovery  factor  is  used  in  calculations. 

1.3.8  Test  Plate  Emissivity 

Radiation  from  the  heated  test  plates  primarily  falls  within 
the  Infrared  range  from  2  to  100  microns.  Plexiglass  has  a  high 
absorbtlvlty  at  these  wavelengths  and  transmits  only  2%  of  the 
Incident  infrared  radiation.  Since  the  test  plates  are  enclosed  by 
the  plexiglass  side  and  top  walls,  a  gray  body  enclosure  radiation 
model  Is  used.  Because  plexiglass  has  a  high  emissivity  of  about 
0.9  and  because  of  the  magnitude  of  the  areas  Involved,  the  general 
gray  body  enclosure  model  simplifies  to  the  special  case  of  a  small 
object  In  a  large  cavity. 
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As  discussed  previously  with  equation  (1.6),  the  radiative 
heat  losses  from  the  test  plates  are  modeled  using 

qr  -  ocA(T4p  -  T4r) 

The  emissivlty  e  of  the  test  plates  Is  very  dependent  on  the  state 
of  oxidation  and  cleanliness  of  the  plates.  Typical  values  of  the 
emissivity  as  quoted  in  various  handbooks  are  -0.05  for  a  polished 
electroplated  nickel  surface  at  23  C,  -0.11  for  an  unpolished 
electroplated  nickel  surface  at  20  C,  and  -0.37  for  heavily  oxi¬ 
dized  nickel  at  200  C.  An  emissivity  typical  of  unpolished 
electroplated  nickel  (0.11)  was  assumed  with  a  bias  limit  of  ± 
0.05. 


1.4  UNCERTAINTY  ANALYSIS  RESULTS 

A  jitter  program  [37]  was  used  to  determine  an  uncertainty 
estimate  for  each  experimentally  determined  Stanton  number.  In  the 
jitter  program,  the  data  reduction  computer  program  was  treated  as 
a  subroutine  and  used  to  approximate  the  partial  derivatives  in 
equation  (1.16)  using  finite  differences.  In  this  numerical  uncer¬ 
tainty  analysis  scheme,  the  first  step  was  to  read  the  original 
experimental  data  from  a  data  file  and  calculate  the  Stanton  number 
for  each  of  the  24  plates  using  the  data  reduction  subroutine. 

Then  the  partial  derivatives  of  the  Stanton  number  with  respect  to 
each  of  the  thirteen  variables  were  determined  numerically  using 
finite  difference  approximations.  The  evaluation  of  a  partial 
derivative  was  accomplished  using  a  loop  through  which  the  original 
value  of  one  variable  was  perturbed  by  a  prespecified  amount  and 
calculating  new  Stanton  numbers  for  all  24  plates  using  the  data 
reduction  subroutine.  The  difference  between  each  new  Stanton 
number  and  the  original  unperturbed  Stanton  number  divided  by  the 
amount  that  the  particular  variable  was  perturbed  resulted  In  a 
value  for  a  partial  derivative.  Before  calculation  of  the  next 
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partial  derivative  the  perturbed  value  of  the  previous  variable  was 
reset  to  the  original  value  and  then  the  next  variable  was  per¬ 
turbed. 

Upon  completion  of  all  partial  derivatives,  equation  (1.16) 
was  evaluated.  Since  precision  limits  In  the  measurements  were 
assumed  negligible,  only  the  bias  limits  for  the  thirteen  variables 
were  required,  and  the  bias  limits  used  were 


Variable 

Bias  limit 

Nominal  Values 

Plate  Temperature 

0.14  C 

45  C 

Rail  Temperature 

0.4  C 

45  C 

Recovery  Temperature 

0.10  C 

30  C 

Wet-Bulb  Temperature 

1.0  C 

27  C 

Barometric  Pressure 

1.0  mm  Hg 

760  mm  Hg 

Recovery  Factor 

0.09 

0.86 

Power 

0.9% 

20-150  W 

Area 

0.03% 

464.5  cm2 

Air  Velocity 

0.4% 

6-70  m/s 

Cpa1r 

0.5% 

1.006  kJ/kg  C 

0.5% 

1.66  kJ/kg  C 

(UATeff 

45% 

0.42  W/C 

Emlsslvlty 

45% 

0.11 

In  addition,  some  of  the  bias  errors  In  temperatures  were 
correlated  since  the  same  calibration  standard  was  used.  The 
correlated  biases  are 


B' 


B' Trail 


0.04  C 


For  the  Stanton  number  data  discussed  In  this  report,  the 
overall  uncertainty  (as  computed  using  equation  (1. 11))  ranged  from 
about  2  percent  to  5  percent,  depending  on  flow  conditions. 
Uncertainty  bars  are  Indicated  on  representative  data  points  when 
the  St  data  are  plotted. 
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Figure  1.2  Power  Measurement  Calibration  Results. 
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ERRTURE  DIFFERENCE!  (C) 


APPENDIX  II 

BOUNDARY  LAYER  PROBE  MEASUREMENTS 


II.l  HOT-WIRE  ANEMOMETRY 

Profiles  of  boundary  layer  mean  velocity  and  turbulence  compo¬ 
nents  were  measured  with  hut-wire  anemometry.  The  hot-wire  instru¬ 
mentation  consisted  of  a  TSI-IFA  100  Intelligent  Flow  Analyzer 
System  with  two  constant  temperature  anemometer  units  and  an 
HP-3437A  high  speed  system  voltmeter  in  the  HP-3054A  Automatic  Data 
Acquisition  and  Control  System  (ADACS)  which  was  linked  to  an 
HP-Series  9000  Model  220  microcomputer.  Two  hot-wire  probes  were 
used:  a  DANTEC  55P05  horizontal,  boundary-layer  type  probe  and  a 
DANTEC  55P02  45°  slant  probe,  both  of  5  micron  diameter  platinum 
plated  tungsten  wire  with  gold  plated  ends.  Each  of  the  two 
anemometer  units  in  the  flow  analyzer  was  dedicated  to  either  the 
horizontal  or  slanted  probe.  Digital  readings  of  the  anemometer 
voltages  were  made  with  the  HP-3437A  high  speed  voltmeter  connected 
directly  to  the  signal  conditioned  anemometer  output.  These  digital 
voltage  readings  were  then  relayed  to  the  microcomputer. 

The  horizontal  hot-wire  probe  and  the  slanted  hot-wire  probe 
are  shown  schematically  with  their  supports  In  Figures  II.l  and 
II. 2,  respectively.  These  probes  are  Identical  to  those  used  by 
Plmenta  1975  [13]  and  Coleman  1976  [14].  Each  probe  was  mounted  in 
a  custom  made  holder  with  a  vertical  micrometer  head  traversing 
mechanism  having  a  resolution  of  0.025  mm.  Both  probe  holders  were 
supported  by  an  xz  traverse  attached  to  a  special  sled  that  spanned 
the  tunnel  and  rested  on  the  side  walls.  Locating  pins  held  the  xz 
traverse  In  place  on  the  special  sled.  Four  hold  down  screws  fixed 
the  sled  to  the  side  walls  of  the  test  section. 

In  use,  both  probes  had  to  be  lowered  to  the  plate  surface  for 
a  reference  height  above  the  plate  to  be  established.  To  prevent 
the  horizontal  wire  from  hitting  the  wall,  a  keel  or  wall  stop  was 
epoxled  to  the  probe  stem  of  the  horizontal  wire.  A  collar  on  the 
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spindle  of  the  slant  wire  probe  provided  a  wall  stop  for  the  slant 
wire.  The  probes  were  lowered  until  electrical  continuity  through 
the  wall  stop  between  the  probe  stem  and  the  plate  surface  was 
established.  Conductive  ink  was  placed  on  the  back  of  the  horizon¬ 
tal  wire  keel  to  provide  electrical  continuity  between  the  keel  and 
probe  stem.  With  the  hot-wire/plate  clearance  provided  by  the  wall 
stops  krowr,  any  desired  vertical  probe  height  above  the  plate 
could  be  set  to  within  0.025  mm.  When  used  with  the  rough  surface, 
a  cylindrical  rod  slightly  longer  than  the  roughness  height  was 
attached  to  the  wall  stops  so  that  the  smooth  portion  of  the  rough 
wall  could  be  used  as  a  reference  position. 

II. 1.1  Horizontal  Hot-Wire  Probe  Measurements 

Measurements  of  the  mean  velocity  and  the  fluctuating  longitu¬ 
dinal  velocity  component  (u‘2)  were  made  with  the  horizontal 
hot-wire.  The  horizontal  wire  was  aligned  with  the  flow  by  match¬ 
ing  etched  marks  on  the  holder  and  horizontal  traversing  mechanism. 
Once  the  probe  was  aligned  with  the  flow  and  freestream  measure¬ 
ments  had  been  taken,  the  boundary  layer  velocity  profile  measure¬ 
ments  began  with  the  probe  keel  starting  Just  above  the  wall  for 
the  first  measurement.  Measurements  were  typically  taken  at  every 
1-2%  of  the  normalized  velocity  (u/U.)  in  the  inner  region  of  the 
boundary  layer  and  every  2-4%  of  the  normalized  velocity  in  the 
outer  region.  At  each  measurement  position,  1000  Instantaneous 
anemometer  output  voltage  readings  were  taken  0.01  seconds  apart 
and  used  to  compute  1000  corresponding  velocities.  A  fourth  order 
least  squares  calibration  equation  was  used  to  convert  anemometer 
voltages  Into  velocities.  The  mean  of  the  1000  computed  velocities 
was  used  as  the  mean  velocity  at  that  location.  The  longitudinal 
velocity  fluctuation  (ul2)  was  taken  as  the  square  of  the  standard 
deviation  (the  variance)  of  the  1000  computed  velocities.  Experi¬ 
ence  showed  that  stable  averages  were  obtained  using  this  number  of 
readings  over  the  10  second  time  period. 
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II. 1.2  Slant  Hot-Wire  Probe  Measurements 

The  slant  wire  was  used  to  determine  the  Reynolds  shear  stress 
factor  (u'v1),  the  normal  velocity  fluctuation  (v'2),  and  the 
transverse  velocity  fluctuation  (wl2).  The  slant  wire  was  mounted 
on  the  rotatable  spindle  of  the  probe  holder  with  Its  prongs  paral¬ 
lel  to  the  mean  flow  direction  at  any  angle  of  rotation.  The 
spinaie  was  rotated  by  a  cable  drive,  which  could  be  operated  with 
the  probe  in  the  tunnel.  A  "lock-drum"  system  with  eight  radially 
drilled  holes  spaced  45°  apart  and  a  spring  loaded  pin  with  a  lever 
located  at  the  top  of  the  vertical  traverse  mechanism  were  used  to 
lock  the  spindle  into  the  desired  orientation  by  fitting  the  pin 
Into  one  of  the  drilled  holes. 

Alignment  of  the  slant  wire  spindle  with  the  mean  flow  was 
done  In  the  freestream.  A  schematic  of  the  slant  wire  geometry  and 
coordinates  Is  given  in  Figure  II. 3.  The  slant  wire  was  placed  in 
the  horizontal  plane  (0-90°, 270°)  and  the  output  of  the  anemometer 
was  noted  for  these  two  probe  orientations.  Alignment  of  the  probe 
spindle  and  slant  wire  with  the  mean  flow  was  adjusted  by  rotating 
the  probe  stem  around  its  y-axls.  The  alignment  of  the  probe  was 
adjusted  back  and  forth  in  small  Increments  across  the  flow  direc¬ 
tion  In  a  iterative  manner  until  the  difference  between  the  elec¬ 
trical  signals  at  probe  rotation  angles  of  0-90°  and  0-270°  was 
3-5mV  from  a  3-5V  signal.  The  corresponding  error  in  the  indicated 
mean  velocity  due  to  misalignment  of  the  slant  wire  probe  was  less 
than  0.12  m/s  for  freestream  velocities  as  high  as  46  m/s. 

To  determine  v'2,  w'2,  and  u'v' ,  the  slant  wire  was  positioned 
approximately  3.3  mm  above  the  surface  of  the  smooth  plate  (or 
about  4  mm  above  the  smooth  surface  of  the  rough  plate)  and  meas¬ 
urements  were  made  at  three  probe  rotation  angles  0-45°,  90°,  and 
135°.  At  each  probe  rotation  angle,  4000  Instantaneous  anemometer 
output  voltage  readings  0.025  seconds  apart  were  taken  and  used  to 
compute  4000  corresponding  effective  velocities  (ufiff).  A  fourth 
order  least  squares  calibration  correlation  was  used  to  convert 
anemometer  voltages  Into  effective  velocities.  The  fluctuating 
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component  of  the  effective  velocity  (u'eff2)  at  each  rotation  angle 
was  taken  as  the  square  of  the  standard  deviation  (the  variance)  of 
the  4000  computed  effective  velocities.  Experience  showed  that 
this  many  readings  taken  over  the  100  second  time  period  provided 
stable  averages.  _ 

The  values  of  u'eff2  at  the  three  slant  wire  probe  orienta¬ 
tions  were  used  in  conjunction  with  the  value  of  u'2  from  the 
horizontal  wire  measurements  at  the  same  y-positlon  to  solve  a 
system  of  three  linear  equations  for  v'2,  w'2,  and  u'v'.  These 
linear  equations  were  generated  with  equation  (II. 1)  evaluated  at 
the  three  probe  rotation  angles  as  discussed  In  detail  by  Coleman 
1976  [14]. 


u'eff 


A  u 


2  +Kv>2 
+  JK  v 


+  4w'2 

4*  w 


+  D 


+  F 


(II. 1) 


The  coefficients  in  this  equation  depend  on  the  orientation  of  the 
probe  with  respect  to  the  flow  coordinates, 

A  -  cos20  +  kj2  sin20 
B  *  (sin20  +  kj2  cos20)  co$20  +  k22  sin20 
C  -  (sin20  +  kj2  cos20)  sin20  +  k22  cos20 
D  -  (1  -  kj2)  sin  20  cos0 
E  -  (si n2^  +  kj2  cos20  -  k22)  sin  20 
F  -  (1  -  kj2)  sin  20  sin0 

where 

0  is  the  probe  rotation  angle  (45°,  90°,  135°) 

0  is  the  wire  slant  angle  (45°) 

The  constants  kj  and  k2,  which  are  known  for  a  DANTEC  55P02  probe, 
depend  on  the  construction  characteristics  of  the  slant  wire  probe 
and  were  taken  to  have  the  values  of 

kj  -  0.2 

k2  -  1.02 

as  previously  done  by  others  (Plmenta  1975  [13]  and  Coleman  1976 
[14],  for  example). 
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II. 1.3  Hot-WI re  Cal Ibratlon 

The  hot-wires  were  calibrated  In  the  test  section  in  the 
uniform  flow  at  the  nozzle  exit.  A  Pitot  probe  placed  2-3  cm  to 
the  side  of  the  hot-wire  and  coupled  to  one  of  the  differential 
pressure  transducers  discussed  in  Appendix  I  was  used  to  determine 
the  freestream  velocity.  The  f reesti  earn  air  velocity  in  the  tunnel 
was  adjusted  to  different  set  points  over  the  desired  measurement 
range  with  the  static  temperature  of  the  freestream  held  constant 
to  within  ±  0.1°C.  At  each  velocity  set  point,  1000  anemometer 
voltage  readings  were  taken  and  averaged.  A  least  squares  analysis 
was  applied  to  the  average  voltages  and  the  velocities  determined 
from  the  Pitot  tube  to  obtain  a  fourth  order  polynomial  calibration 
equation  for  each  wire. 

In  practice,  the  resistance  of  the  heated  wire  during  measure¬ 
ments  (the  operating  resistance)  is  set  above  the  resistance  of  the 
unheated  wire  (the  cold  probe  resistance)  to  satisfy  the  desired 
overheat  ratio  requirements.  The  cold  probe  resistance  is  tempera¬ 
ture  dependent  and  must  be  remeasured  and  re-entered  into  the 
anemometer  unit  if  the  freestream  temperature  of  the  test  air 
changes  more  than  about  ±  0.2°C  from  the  conditions  of  calibration. 
The  required  operating  resistance  used  was  always  3  ohms  higher 
than  the  cold  probe  resistance  at  all  operating  conditions.  A 
comparison  of  calibration  equations  from  a  calibration  done  at  a 
freestream  temperature  of  25.6°C  (78° F)  and  a  calibration  done  at  a 
freestream  temperature  of  29.9°C  (84°F)  showed  that  the  velocities 
computed  with  the  two  equations  differed  by  less  than  1%  over  the 
entire  calibration  range.  Since  these  calibrations  had  different 
cold  probe  resistances  and  the  operating  resistances  were  set  by 
adding  3  ohms  to  the  cold  probe  resistance.  It  was  concluded  that 
calibrations  of  the  hot-wires  at  one  freestream  temperature  would 
be  valid  over  a  small  temperature  range  (±  3°C)  as  long  as  the 
difference  in  the  probe  resistances  was  held  at  3  ohms  and  the 


freestream  temperature  of  the  test  air  was  not  allowed  to  vary  more 
than  about  ±  0.1°C  while  data  were  taken  for  a  profile  at  a  par¬ 
ticular  x-locatlon. 

II. 1.4  Uncertainties 

Uncertainties  in  the  hot-wire  measurements  were  estimated  by 
considering  the  uncertainties  in  the  Pitot-determined  velocity,  in 
curvefitting  of  the  calibration  data,  In  anemometer  adjustments, 
and  in  probe  alignment;  the  observed  scatter  in  the  data;  and  the 
suggested  uncertainties  given  In  the  literature  (Kline,  Cantwell, 
and  LI  1  ley  1981  [56],  Yavuzkurt  1982  [57]).  The  order  of  the 
overall  uncertainties  (bias  and  precision)  associated  with  the 
hot-wire  measurements  are:  u,  2%;  u'2,  5%;  v'2,  15%;  w'2,  10%;  and 
uV,  10%. 

1 1. 2  THERMOCOUPLE  PROBE 

Time  mean  temperatures  In  the  boundary  layer  were  measured 
using  a  Type  E  (chromel-constantan)  butt-welded  thermocouple  probe 
similar  in  design  to  that  of  Blackwell  1972  [58],  The  output  of 
the  thermocouple  is  in  the  millivolt  range  and  is  measured  by  the 
ADACS. 

The  thermocouple  probe  holder,  which  is  almost  identical  to 
the  horizontal  wire  holder,  was  mounted  on  the  same  support  sled 
used  with  the  hot-wires.  A  schematic  of  the  thermocouple  probe  and 
its  support  is  shown  in  Figure  II. 4. 

The  thermocouple  calibration  was  done  in  a  Blue  M  Model 
MR-3210A-1  constant  temperature  bath,  and  the  Hewlett-Packard 
quartz  thermometer  described  In  Appendix  I  was  utilized  to  measure 
the  bath  temperature.  The  calibration  water  bath  was  In  continuous 
movement  due  to  an  automatic  stirrer,  and  the  risk  of  breaking  the 
fine  thermocouple  wire  was  large  if  the  probe  was  placed  directly 
into  it.  Besides,  water  could  deposit  some  residue  on  the  wire 
surface  and  the  prongs,  which  could  Influence  the  thermocouple 
temperature  response.  Moreover,  the  water  temperature  close  to  the 
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thermocouple  could  not  be  accurately  monitored  by  the  quartz  ther¬ 
mometer.  To  alleviate  those  difficulties  It  was  decided  to  Insert 
the  wire  Into  a  jar  filled  with  alcohol,  which  was  placed  In  the 
water  bath.  The  quartz  thermometer  was  also  positioned  In  the  jar 
next  to  the  thermocouple  probe  so  that  It  would  encounter  the  same 
conditions.  To  prevent  any  air  current  from  convecting  heat  to  or 
from  the  alcohol  surface,  the  opening  of  the  jar  was  sealed.  The 
time  constant  of  the  jar  was  also  accounted  for  by  waiting  an  hour 
after  the  water  bath  temperature  had  reached  the  steady  state 
condition  before  proceeding  with  the  temperature  measurement. 

Calibration  was  performed  for  temperatures  between  23°C-39°C 
using  four  points  over  this  range.  The  thermocouple  probe  voltage 
outputs,  as  measured  by  the  AOACS,  were  converted  to  temperatures 
using  the  HP  system  software  package.  The  temperature  of  the 
reference  junction  (the  isothermal  terminal  block)  required  for 
software  compensation  Is  established  by  the  ADACS  via  a  temperature 
transducer  which  provides  a  100  mV/°C  output  voltage.  This  soft¬ 
ware  performs  the  voltage-temperature  conversion  by  dividing  the 
thermocouple  characteristic  curve  Into  eight  sectors  and  approxi¬ 
mating  each  sector  by  a  third  order  nested  polynomial.  The  tem¬ 
peratures  measured  by  the  thermocouple  were  compared  with  the 
temperatures  from  the  quartz  thermometer  over  the  above  mentioned 
calibration  range.  The  departure  of  thermocouple  temperatures, 
using  the  software  package  for  conversion,  from  the  corresponding 
temperatures  obtained  by  the  quartz  thermometer  was  less  than  ± 
0.08°C. 
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APPENDIX  III 

SKIN  FRICTION  COEFFICIENT  DETERMINATION 


Skin  friction  coefficients  in  this  study  were  determined  in 
two  ways:  (1)  from  hot-wire  measurements  of  Reynolds  shear  stress 
and  mean  velocity  profiles,  and  (2)  from  Preston  tube  measurements. 
These  skin  friction  determination  methods  are  discussed  briefly 
below  with  detailed  references  given  for  further  Information. 

III.l  HOT-WIRE  METHOD 

The  first  means  of  determining  the  skin  friction  coefficients 
in  this  study  was  the  hot-wire  data  method.  (This  procedure  Is  the 
only  method  used  for  skin  friction  determination  In  the  rough  wall 
boundary  layers.)  Skin  friction  coefficients  were  determined  using 


Cf  v  3u 


Yl 


+  ij/lfc)  dy^  (III*1) 

To  derive  equation  (III.l),  the  Reynolds-averaged  momentum 
equation  (incorporating  the  usual  boundary  layer  assumptions)  and 
the  continuity  equation  are  integrated  from  the  plate  surface  to  a 
position  Yj  in  the  boundary  layer.  Calculations  of  Cf  are  made 
using  equation  (III.l),  mean  velocity  profiles  measured  over  adja¬ 
cent  test  plates  (AX  -  10.2  cm),  and  the  measured  value  u' v'  at 
y-Yj.  The  position  Yj  was  about  3.3  mm  for  smooth  wall  studies  and 
4  mm  for  rough  wall  studies.  Further  details  of  hot-wire  data  de¬ 
terminations  of  the  skin  friction  coefficient  are  given  by  Andersen 
1972  [59],  Pimenta  1975  [13],  and  Coleman  1976  [14]. 
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The  estimated  uncertainty  in  the  hot-wire  determined  skin 
friction  coefficients  In  this  study  Is  about  ±10-12%.  This  esti¬ 
mate  Is  based  on  smooth  wall  comparisons  of  the  hot-wire  determined 
values  of  Cf  with  accepted  correlations  and  with  the  values 
determined  using  the  Preston  tube. 

I I I. 2  PRESTON  TUBE  METHOD 

The  method  of  Preston  1954  [60]  for  determining  the  skin  fric¬ 
tion  in  turbulent  boundary  layer  flows  uses  a  simple  Pitot  tube 
(Preston  tube)  resting  on  the  surface  and  depends  upon  the  assump¬ 
tion  of  a  universal  inner  law  (law  of  the  wall)  common  to  smooth 
wall  boundary  layer  flows.  The  difference  between  the  total  pres¬ 
sure  at  the  Preston  tube  of  1.6  mm  (1/16  in)  inside  diameter  and 
the  undisturbed  static  pressure  at  a  pressure  tap  in  the  test 
section  sidewall  at  the  same  x-location  was  measured  with  the 
pressure  transducers  described  in  Appendix  I.  This  difference  in 
pressure  was  then  used  in  conjunction  with  the  calibration  equa¬ 
tions  as  given  by  Patel  1965  [61]  to  solve  for  the  local  skin  fric¬ 
tion  coefficient. 

The  Preston  tube  method  of  determining  skin  friction  coeffi¬ 
cients  is  only  of  use  for  the  smooth  plate  case  because  calibra¬ 
tions  have  only  been  made  for  smooth  walls.  However,  skin 
friction  coefficient  determinations  made  by  Preston's  method  are 
quick  and  easy  to  make.  The  relatively  low  uncertainty  (less  than 
about  6%)  in  the  skin  friction  coefficients  determined  with 
Preston's  method  made  this  method  very  important  in  establishing 
the  qualification  of  the  test  facility  and  the  correctness  of  the 
hot-wire  based  skin  friction  determination  method. 
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